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Abstract 

Some aspects of physics on interacting fermions in ID are discussed in a tutorial- 
like manner. We begin by showing that the non-analytic corrections to the Fermi- 
liquid forms of thermodynamic quantities result from essentially ID collisions em- 
bedded into a higher-dimensional phase space. The role of these collisions in- 
creases progressively as dimensionality is reduced until, finally, they lead to a 
breakdown of the Fermi liquid in ID. An exact solution of the Tomonaga-Luttinger 
model, based on the Ward identities, is reviewed in the fermionic language. Tunnel- 
ing in a ID interacting systems is discussed first in terms of the scattering theory for 
interacting fermions and then via bosonization. Universality of conductance quan- 
tization in disorder-free quantum wires is discussed along with the breakdown of 
this universality in the spin-incoherent case. A difference between charge (univer- 
sal) and thermal (non-universal) conductances is explained in terms of Fabry-Perrot 
resonances of charge plasmons. 



Fundamental aspects of electron correlations and quantum transport 5 
1. Introduction 

The theory of interacting fermions in one dimension (ID) has survived several 
metamorphoses. From what seemed to be a purely mathematical exercise up until 
the 60s, it had evolved into a practical tool for predicting and describing phenom- 
ena in conducting polymers and organic compounds- which were the ID systems 
of the 70s. Beginning from the early 90s, when the progress in nanofabrication 
led to creation of artificial ID structures-quantum wires and carbon nanotubes, 
the theory of ID systems started its expansion into the domain of mesoscopics; 
this trend promises to continue in the future. Given that there is already quite a 
few excellent reviews and books on the subject [1]-[10] , I should probably be- 
gin with an explanation as to what makes this review different from the others. 
First of all, it is not a review but-being almost a verbatim transcript of the lec- 
tures given at the 2004 Summer School in Les Houches-rather a tutorial on some 
(and definitely not all) aspects of ID physics. A typical review on the subject 
starts with describing the Fermi Liquid (FL) in higher dimensions with an aim of 
emphasizing the differences between the FL and its ID counter-part -Luttinger 
Liquid (LL). My goal-if defined only after the manuscript was written-was rather 
to highlight the similarities between higher-D and ID systems. The progress in 
understanding of ID systems has been faciUtated tremendously and advanced to 
a greater detail, as compared to higher dimensions, by the availability of exact or 
asymptotically exact methods (Bethe Ansatz, bosonization, conformal field the- 
ory), which typically do not work too well above ID. The downside part of this 
progress is that ID effects, being studied by specifically ID methods, look some- 
what special and not obviously related to higher dimensions. Actually, this is 
not true. Many effects that are viewed as hallmarks of ID physics, e.g., the sup- 
pression of the tunneling conductance by the electron-electron interaction and 
the infrared catastrophe, do have higher-D counter-parts and stem from essen- 
tially the same physics. For example, scattering at Friedel oscillations caused by 
tunneling barriers and impurities is responsible for zero-bias tunneling anoma- 
lies in all dimensions 16]. The difference is in the magnitude of the effect 
but not in its quahtative nature. Following the tradition, I also start with the FL 
in Sec. 2, but the main message of this Section is that the difference between 
D = 1 and Z? > 1 is not all that dramatic. In particular, it is shown that the 
well-known non-analytic corrections to the FL forms of thermodynamic quan- 
tities (such as a venerable In T-correction to the linear-in-T specific heat in 
3D) stem from rare events of essentially ID collisions embedded into a higher- 
dimensional phase space. In this approach, the difference between D = 1 and 
£> > 1 is quantitative rather than quahtative: as the dimensionality goes down, 
the phase space has difficulties suppressing the small-angle and 2/2^?— scattering 
events, which are responsible for non-analyticities. The special point when these 



6 



Dmitrii L. Maslov 



events go out of control and start to dominate the physics happens to be in ID. 

This theme is continued in Sec. 5, where scattering from a single impurity em- 
bedded into a ID system is analyzed in the fermionic language, following the 
work by Yue, Matveev, Glazman [11]. The drawback of this approach-the per- 
turbative treatment of the interaction-is compensated by the clarity of underlying 
physics. Another feature which makes these notes different from the rest of the 
Uterature in the field is that the description goes in terms of the original fermions 
for quite a while (Sees. 2 through 5) , whereas the weapon of choice of all ID 
studies-bosonization- is invoked only at a later stage (Sec. 6 and beyond). The 
rationale-again, a post-factum one-is two-fold. First, ID systems in a meso- 
scopic environment-which are the main real-life application discussed here- are 
invariably coupled to the outside world via leads, gates, etc. As the outside world 
is inhabited by real fermions, it is sometimes easier to think of, e.g., both the 
interior and exterior a quantum wire coupled to reservoirs in terms of the same 
elementary quasi-particles. Second, after 40 years or so of bosonization, what 
could have been studied within a model of fermions with linearized dispersion 
and not too strong interaction-and this is when bosonization works-was proba- 
bly studied. (As all statements of this kind, this is one is also an exaggeration.) 
The last couple of years are characterized by a growing interest in either the ef- 
fects that do not occur in a model with linearized dispersion, e.g., Coulomb drag 
due to small-momentum transfers [17| and energy relaxation, or situations when 
strong Coulomb repulsion does not permit linearization of the spectrum at any 
energies [19, 20, 21]. Experiment seems to indicate that the Coulomb repulsion 
is strong in most systems of interest, thus studies of a strongly-coupled regime 
are quite timely. Once the assumption of the linear spectrum is abandoned, the 
beauty of a bosonized description is by and large lost, and one might as well 
come back to original fermions. Sec. 6 is devoted to transport in quantum wires, 
mostly in the absence of impurities. The universality of conductance quantiza- 
tion is explained in some detail, and is followed by a brief discussion of the recent 
result by Matveev [19], who showed that incoherence in the spin sector leads to 
a breakdown of the universality at higher temperatures (Sec. 7.4). Also, a differ- 
ence in charge (universal) and thermal (non-universal) transport-emphasized by 
Fazio, Hekking, and Khmelnitskii [22]- in addressed in Sec.7.5. What is missing 
is a discussion of transport in a disordered (as opposed to a single-impurity) case. 
However, this canonically difficult subject, which involves an interplay between 
locaUzation and interaction, is perhaps not ready for a tutorial-Uke discussion at 
the moment. (For a recent development on this subject, see Ref.[18].) 

Even a brief inspection of these notes shows that the choice between making 
them comprehensive or self-contained was made for the latter with a focus on a 
relatively small number of topics. It is quite easy to see what is missing: there 
is no discussion of lattice effects, bosonization is introduced without the Klein 
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factors, the sine-Gordon model is not treated in depth, chiral Luttinger liquids 
are not discussed at all, and the list goes on. The discussion of the experiment 
is scarce and perfunctory. However, the few subjects that are discussed are pro- 
vided with quite a detailed-perhaps somewhat excessively detailed- treatment, 
so that a reader may not feel a need to consult the reference list too often. For 
the same reason, the notes also cover such canonical procedures as the pertur- 
bative renormalization group in the fermionic language (Sec. 4) and elementary 
bosonization (Sec. 6), which are discussed in many other sources and a reader 
already familiar with the subject is encouraged to skip them. 

Also, a relatively small number of references (about one per page on average) 
indicates once again that this is not a review. The choice of cited papers is sub- 
jective and the reference list in no way pretends to represent a comprehensive 
bibliography to the field. My apologies in advance to those whose contributions 
to the field I have failed to acknowledge here. 

/I = fcs = 1 through out the notes, unless specified otherwise. 



2. Non-Fermi liquid features of Fermi liquids: ID physics in higher dimen- 
sions 



One often hears the statement that, by and large, a Fermi liquid (FL) is just a 
Fermi gas of weakly interacting quasi-particles; the only difference being the 
renormalization of the essential parameters (effective mass, g— factor) by the in- 
teractions. What is missing here is that the similarity between the FL and Fermi 
gas holds only for leading terms in the expansion of the thermodynamic quan- 
tities (specific heat C(T), spin susceptibiUty Xs> etc.) in the energy (tempera- 
ture) or spatial (momentum) scales. Next-to-leading terms, although subleading, 
are singular (non-analytic) and, upon a closer inspection, reveal a rich physics 
of essentially ID scattering processes, embedded into a high-dimensional phase 
space. 

In this chapter, 1 will discuss the difference between "normal" processes which 
lead to the leading, FL forms of thermodynamic quantities and "rare" ID pro- 
cesses which are responsible for the non-analytic behavior. We will see that the 
role of these rare processes increases as the dimensionality is reduced and, even- 
tually, the rare processes become the norm in ID, where the FL breaks down. 

In a Fermi gas, thermodynamic quantities form regular, analytic series as func- 
tion of either temperature, T, or the inverse spatial scale (bosonic momentum q) 
of an inhomogeneous magnetic field. For T <C Ep, where Ep is the Fermi 
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energy, and q <^ kp, where fci? is the Fermi momentum, we have 



C{T)/T 
Xs {T,q = 0) 
Xs {T = 0,q) 



-f + aT'^ + bT^ + ...; 
X'l{0) + eq^ + fq^ + ..., 



(2.1a) 
(2.1b) 
(2.1c) 



where 7 is the Sommerfeld constant, x° is the static, zero-temperature spin sus- 
ceptibility (which is finite in the Fermi gas), and a . . . / are some constants. Even 
powers of T occur because of the approximate particle-hole symmetry of the 
Fermi function around the Fermi energy and even powers of q arise because of 
the analyticity requirement ' 

Our knowledge of the interacting systems comes from two sources. For a sys- 
tem with repulsive interactions, one can assume that as long as the strength of the 
interaction does not exceed some critical value, none of the symmetries (transla- 
tional invariance, time-reversal, spin-rotation, etc.), inherent to the original Fermi 
gas, are broken. In this range, the FL theory is supposed to work. However, the 
FL theory is an asymptotically low-energy theory by construction, and it is really 
suitable only for extracting the leading terms, corresponding to the first terms 
in the Fermi-gas expressions (2.1a-2.1c). Indeed, the free energy of a FL as an 
ensemble of quasi-particles interacting in a pair-wise manner can be written as 



where Fq is the ground state energy, 5nk is the deviation of the fermion occu- 
pation number from its ground-state value, and fk,k' is the Landau interaction 
function. As Suk is of the order of T/Ep, the free energy is at most quadratic 
in T, and therefore the corresponding C{T) is at most hnear in T. Consequently, 
the FL theory-at least, in the conventional formulation-claims only that 



The expressions presented above are valid in all dimensions, except for D = 2 with quadratic 
dispersion. There, because the density of states (DoS) does not depend on energy, the leading correc- 
tion to the 7T— term in C(T) is exponential in Ep /T and Xs does not depend on q for q < 2kp. 
However, this anomaly is removed as soon as we take into account a finite bandwidth of the electron 
spectrum, upon which the universal (T^" and g^" ) behavior of the series is restored. 



[25] 




(2.2) 



C*(T)/T = 7*; 
X*AT,q) = x:(0) 
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Fig. 1. Combined "diagram of knowledge", x-axis: energy scale (given by temperature T, bosonic 
momentum Q, magnetic field H) in appropriate units, y-axis: interaction strength. Fermi liquid 
works for not necessarily weak interactions (but smaller than the critical value for an instability of 
the ground state denoted by the red dot) but at the lowest energy scales. Microscopic models work 
for weak interactions but for arbitrary energies. 

where 7* and xl (0) differ from the corresponding Fermi-gas values, and does 
not say anything about higher-order terms ^. 

Higher-order terms in T or g can be obtained within microscopic models 
which specify particular interaction and, if an exact solution is impossible-which 
is always the case in higher dimensions- employ some kind of a perturbation the- 
ory. Such an approach is complementary to the FL: the former nominally works 
for weak interactions ^ but at arbitrary temperatures, whereas FL works both 
for weak and strong interactions, up to some critical value corresponding to an 
instability of some kind, e.g., a ferromagnetic transition, but only in the low- 
temperature limit. In the { temperature, interaction} plane, the validity regions 
of these two approaches are two strips running along the two axes (cf. Fig. 1). 
For weak interactions and at low temperatures, the regions should overlap. 

Microscopic models (Fermi gas with weak repulsion. Coulomb gas in the 
high-density limit, electron-phonon interaction, paramagnon model, etc.) show 
that the higher-order terms in the specific heat and spin susceptibility are non- 

^ Strictly speaking, non-analytic terms in C(T) can be obtained from the free energy (2.2) by 
taking into account the non-analytic terms in the quasi-particle spectrum, see Ref. [29]b. 

Some results of the perturbation theory can be rigorously extended to an infinite order in the 
interaction, and most of them can be guaranteed to hold even if the interactions are not weak. 
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analytic functions of T and g [26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38]. 
For example. 



C{T)/T = 73 -agT^ In T (3D); (2.3a) 

C{T)/T = 72-a2T(2D); (2.3b) 

Xs{q) = x.(0) + /33g'lng-i(3D); (2.3c) 

Xsid) = X«(0)+/32|g| (2D), (2.3d) 



where all coefficients are positive for the case of repulsive electron-electron in- 
teraction 

As seen from Eqs. (2.3a-2.3d), the non-analyticities become stronger as the 
dimensionality is reduced. The strongest non-analyticity occurs in ID, where- 
as far as single-particle properties are concerned-the FL breaks down: 

C{T)/T = 71 +ailnT (ID); 
Xil) = Xo + /3iln|«| (ID). 

It turns out that the evolution of the non-analytic behavior with the dimen- 
sionality reflects an increasing role of special, almost ID scattering processes in 
higher dimensions. Thus non-analyticities in higher dimensions can be viewed 
as precursors of ID physics for D > 1. 

It is easier to start with the non-analytic behavior of a single-particle prop- 
erty, the self-energy, which can be related to the thermodynamic quantities via 
standard means [23] (see also appendix A). Within the Fermi Uquid, 

ReS^(£,fc) = -Ae + Bik + ■■■ (2.4a) 
-ImS^(£,fc) = C(e2 -H TT^T^) -h . . . (2.4b) 

Expressions (2.4a) and (2.4b) are equivalent to two statements: i) quasi-particles 
have a finite effective mass near the Fermi level 

''Notice that not only the functional forms but also the sign of the q— dependent term in the 
spin susceptibility is different for free and interacting systems. "Wrong" sign of the q— dependent 
corrections has far-reaching consequences for quantum critical phenomena. For example, it precludes 
a possibility of a second-order, homogeneous quantum ferromagnetic phase transition in an itinerant 
system [39]. What is possible is either a first-order transition or ordering at finite q, e.g. heUcal 
structure. In ID, a homogeneous ferromagnetic state is forbidden by the Lieb-Mattis theorem [46], 
which states that the ground state of ID fermions, interacting via spin-independent but otherwise 
arbitrary forces, is non-magnetic. One could speculate that the non-analyticities in higher dimensions 
indicate the existence of a higher-D version of the Lieb-Mattis theorem. Certainly, this does not 
mean that ferromagnetism does not exist in higher dimensions (it is hard to deny the existence of, 
e.g., iron). However, ferromagnetism may not exist in models dealing only with itinerant electrons in 
continuum. 
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Fig. 2. Self-energy to first order in tlie interaction with a dynamic bosonic field. 



m* = mo 



A + 1 
B + V 



and ii) damping of quasiparticles is weak: the level width is much smaller than 
the typical quasi-particle energy 



Landau's argument for the (or T^) behavior of ImS^ relies on the Fermi 
statistics of quasiparticles and on the assumption that the effective interaction is 
screened at large distances [23]. It requires two conditions. One condition is ob- 
vious: the temperature has to be much smaller than the degeneracy temperature 
Tp = kpVp, where Vp is the renormahzed Fermi velocity. The other condition 
is less obvious: it requires inter-particle scattering to be dominated by processes 
with large (generically, of order kp) momentum transfers. Once these two con- 
ditions are satisfied, the self-energy assumes a universal form, Eqs. (2.4a) and ( 
2.4b), regardless of a specific type of the interaction (e-e, e-ph) and dimensional- 
ity. To see this, let's have a look at ImS^ (e, k) due to the interaction with some 
"boson" (Fig. 2). 

The wavy fine in Fig.2 can be, e.g., a dynamic Coulomb interaction, phonon 
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propagator, etc. On the mass shell {e = ^fc) at T = and for £ > 0, we have ^ 

ImS^ (e) = ?g:rr j I d'^qlmG^ (e - w, k - q) ImV^ {oj, q) . 

(27r) Jo J 

(2.5) 

The constraint on energy transfers (0 < w < e) is a direct manifestation of 
the Pauli principle which limits the number of accessible energy levels. In real 
space and time, V (r, t) is a propagator of some field which has a classical limit 
(when the occupation numbers of all modes are large). Therefore, V (r, t) is a 
real function, hence ImV is an odd function of w. I wiU make this fact expUcit 
writing ImV as 

ImV^ {uj,q) =ojW {\uj\,q) . 

Now, suppose that we integrate over q and the result does not depend on co. Then 
we innmediately get 

-ImE^ (e) ~ C / (Ljcj ~ C£^ 
Jo 

where C is the result of the q— integration which contains all the information 
about the interaction. Once we got the e^-form for ImS^ (e) , the e- term in 
ReS^ (e) follows immediately from the Kramers-Kronig transformation, and we 
have a Fermi-liquid form of the self-energy regardless of a particular interaction 
and dimensionality. Thus a sufficient condition for the Fermi liquid is the sepa- 
rability of the frequency and momentum integrations, which can only happen if 
the energy and momentum transfers are decoupled. 

Now, what is the condition for separability? As a function of q, W has at 
least two characteristic scales. One is provided by the internal structure of the 
interaction (screening wavevector for the Coulomb potential, Debye wavevector 
for electron-phonon interaction, etc.) or by kp, whichever is smaller. This scale 
(let's call it Q) does not depend on co. Moreover, as |a;| is bounded from above 
by £, and we are interested in the limit £ — > 0, one can safely assume that Q ^ 
\uj \ /vp. The role of Q is just to guarantee the convergence of the momentum 
integral in the ultraviolet, that is, to ensure that for q ^ Q the integrand falls 
off rapidly enough. Any physical interaction will have this property as larger 
momentum transfer will have smaller weight. The other scale is \ /vp- Now, 

^To get Eq. (2.5), one can start with the Matsubara form of diagram Fig. 2, convert the Matsubara 
sums into the contour integrals, use the dispersion relation 

which is valid for any retarded function, and take the limit T — > 0. 
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let's summarize this by re-writing \mV in the following scaling form 



where C/ is a dimensionless function and the factor was singled out to keep 
the units right. 

In the perturbation theory, the Green's function in (2.5) is a free one. Assum- 
ing the free-electron spectrum ^fc = (fc^ — kp)/2m, 

ImG^ ^£ — LO,k — ^ = —ttS [s — to — +Vk ■ q — /2rn) . 

On the mass shell, 

ImG^ (^s - u>,k - qj \e=^^^ = -ttS (u> -Vk ■ q + q^ /2m) . 

The argument of the delta-function simply expresses the energy and momentum 
conservation for a process e ^ e — uj,k ^ k — q. The angular integral involves 
only the delta-function. For any D, this integral gives 

(()(... ))n = —Ad 



Vpq \ vpq 

where Vk was replaced by ii^^ because all the action takes place near the Fermi 
surface. For = 3 and D = 2, 

Asix) = 20(1 -\x\); 

2e{i-\x\) 



A2ix) = 



The constraint on the argument oi Ad is purely geometric: the magnitude of the 
cosine of the angle between k and ^has to be less then one. For power-counting 
purposes, function Ad has a dimensionality of 1. Therefore, its only role is to 
provide a lower cut-off for the momentum integral. Then, by power counting 

ImS« (e) ^ r cLouj I dqq^-^U (^,^). (2.6) 

vfQ'^ Jo Jq>\u,\/vp \'"Fq QJ 

Now, if the integral over q is dominated by g ~ Q and is convergent in the in- 
frared, one can put = in this integral. After this step, the integrals over oj and 
q decouple. The w— integral gives £^ regardless of the nature of the interaction 
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Fig. 3. a) and b) Non-trivial second order diagrams for tlie self-energy, c) Same diagrams as in a) 
and b) re-drawn as a single "sunrise" diagram, d) Diagrams relevant for non-analytic terms in the 
self-energy, e) Kinematics of scattering in a polarization bubble: the dynamic part 11 oc uj/v^q 
comes from the processes in which the internal fermionic momentum (p) is almost perpendicular to 
the external bosonic one (q). 
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and dimensionality whereas the q— integral supplies a prefactor which entails all 
the details of the interaction 



ImS^ (e) = Cd 



VfQ' 



For example, for a screened Coulomb interaction in the weak-coupling (high- 
density) limit Q = K, where k is the screening wavevector, we have in 3D 

2 2 

-ImS^ (e) = — — — . 

Now we can formulate a sufficient (but not necessary) condition for the Fermi- 
liquid behavior. It will occur whenever if kinematics of scattering is such that 
the typical momentum transfers are determined by some internal and, what is 
crucial, oj— independent scale, whereas the energy transfers are of order of the 
quasi-particle energy (or temperature). Excluding special situations, such as 
the high-density hmit of the Coulomb interaction, Q is generically of order of 
the ultraviolet range of the problem ~ kp-^n other words, isotropic scattering 
guarantees a e^- behavior. Small-angle scattering with typical angles of order 
e/vp <^ Q <^ kp gives this behavior as well. 

The s^— result seems to be quite general under the assumptions made. When 
and why these assumptions are violated? 

A long-range interaction, associated with small-angle scattering, is known to 
destroy the FL. For example, transverse long-range (current-current [44J or gauge 
[45]) interactions, which-unlike the Coulomb one-are not screened by electrons, 
lead to the breakdown of the Fermi liquid. However, the current-current inter- 
action is of the relativistic origin and hence does the trick only at relativistically 
small energy scales, whereas the gauge interaction occurs only under special cir- 
cumstances, such as near half-filhng or for composite fermions. What about a 
generic case when nothing of this kind happens? It turns out that even if the 
bare interaction is of the most benign form, e.g., a delta-function in real space, 
there are deviations from a (perceived) FL behavior. These deviations get amph- 
fied as the system dimensionality is lowered, and, eventually, lead to a complete 
breakdown of the FL in ID. 

A formal reason for the deviation from the FL-behavior is that the argument 
which led us to the e^-term is good only in the leading order in to/qvp. Recall 
that the angular integration gives us q~^ factors in all dimensions, and, to arrive 
at the £^ result we put a; = in functions Ajj and U. If we want to get a next 
term in e, then we need to expand U and Ainu). Had such expansions generated 
regular series, ImS^ would have also formed regular series in e^: ImS^ = 
ae^ + + ce^ + However, each factor of ui comes with q~^, so that no 
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matter how high the dimensionahty is, at some order of uj/ypq we are bound to 
have an infrared divergence. 

2.1. Long-range effective interaction 

Let's look at the simplest case of a point-like interaction. A frequency depen- 
dence of the self-energy arises already at the second order. At this order, two 
diagrams in Fig. 3 are of interest to us. For a contact interaction, diagram b) is 
just -1/2 of a) (which can be seen by integrating over the fermionic momentum p 
first), so we will lump them together. Two given fermions interact via polarizing 
the medium consisting of other fermions. Hence, the effective interaction at the 
second order is just proportional to the polarization bubble 

Let's focus on small angle-scattering first: q <C 2kp. It turns out that in all 
three dimensions, the bubble has a similar form (see Appendix Appendix A for 
an exphcit derivation of this result) 

-Imn«(a;,5) = vd—Bd ( — ) , (2.7) 
Vpq \vpqj 

where vd = aDmkp~^ is the DoS in D dimensions [03 = (27r)^^, 02 = 
(27r)^^, ai = 1/2tt] and Bd is a dimensionless function, whose main role is 
to impose a constraint iv < Vpq in 2D and 3D and w = Vpq in ID. Eq. (2.7) 
entails the physics of Landau damping. The constraint arises because collective 
excitations-charge- and spin-density waves- decay into particIe-hoIe pairs. De- 
cay occurs only if bosonic momentum and frequency (g and to) are within the 
particle-hole continuum (cf. Fig. 4). For D > 1, the boundary of the continuum 
for small u) and q\suj = vpq, hence the decay takes place \f oj < vpq. The 
rest of Eq. (2.7) can be understood by dimensional analysis. Indeed, 11^ is the 
retarded density-density correlation function; hence, by the same argument we 
applied to ImT^^, its imaginary part must be odd in uj. For q <C kp, the only 
combination of units of frequency is vpq, and the frequency enters as ui/vpq. 
Finally, a factor ud makes the overall units right. In ID, the difference is that the 
continuum shrinks to a single line lu = vpq, hence decay of collective excitations 
is possible only on this line. In 3D, function B3 is simply a 9— function 

Imn^(a;,g) = -u3—e{q-\ij/vF\). 
Vpq 

Next-to-leading term in the expansion of ImS^ in e comes from retaining the 
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(a) 



(b) 



Fig. 4. Particle-hole continua for D > I (left) and D = 1 (right). For the ID case, only half of the 
continuum {q > 0) is shown. 




Ill /v„ 



q~ loil /Vp 



Iq - 2kpl~ |ai| /Vp 





Fig. 5. Kinematics of scattering, a) "Any-angle" scattering. Momentum transfer q is of order of the 
intrinsic scale of the interaction or kp, whichever is smaller, and is independent of the energy transfer, 
uj, which is of order of the initial energy e. This process contributes regular terms to the self-energy, 
b) Dynamic forward scattering: q ~ \Lo\/vp. c) Dynamic backscattering: \q — 2kp\ ~ \u)\/vp. 
Processes b) and c) are responsible for the non-analytic terms in the self-energy. 
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lower limit in the momentum integral of Eq. (2.6), upon which we get 

-ImS^ ~ U'^mkp / dw / dqq^ 

Jo Jio/vp vpq vpq 

^ — 5— / dcoui 
K Jo 

~ ae^ — b \e\^ . 

The first term in the square brackets is the FL contribution that comes from 
q Q. The second term is a correction to the FL coming from q ^ cj/vp- 
Thus, contrary to a naive expectation an expansion in e is non-analytic. The frac- 
tion of phase space for small-angle scattering is small-most of the self-energy 
comes from large-angle scattering events (g ~ Q); but we already start to see the 
importance of the small-angle processes. Applying Kramers-Kronig transforma- 
tion to the non-analytic part (|£|^) in ImS^, we get a corresponding non-analytic 
contribution to the real part as 

(ReS^) cx In \s\ . 

V /non— an ' ' 

Correspondingly, specific heat which, by power counting, is obtained from ReS^ 
by replacing each e by T, also acquires a non-analytic term^ 

C(T) =73T + /33r'Mnr. 

This is the familiar T"* InT term, observed both in He^ [40] and metals [41] 
(mostly, heavy-fermion materials) ^. 

In 2D, the situation is more dramatic. The q— integral diverges now logarith- 

*One has to be careful with the argument, as a general relation between C(T) and the single- 
particle Green's function [23] involves the self-energy on the mass shell. In 3D, the contribution 
to S from forward scattering, as defined in Fig. 6, vanishes on the mass shell; hence there is no 
contribution to C{T) [50]. The non-analytic part of C(T) is related to the backscattering part of the 
self-energy (scattering of fermions with small total momentum), which remains finite on the mass 
shell. That forward scattering does not contribute to non-analyticities in thermodynamics is a general 
property of all dimensions, which can be understood on the basis of gauge-invariance [42] . 

'The T"* In T-term in the specific heat coming from the electron-electron interactions is often re- 
ferred to in the literature as to the "spin-fluctuation" or "paramagnon" contribution [27, 28]. Whereas 
it is true that this term is enhanced in the vicinity of a ferromagnetic (Stoner) instability, it exists even 
far way from any critical point and arises already at the second order in the interaction [29]. 



U) 

kp — — 

FL ^-v-' 

beyond FL 
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mically in the infrared: 



-ImE-^ 



(a;) 




'^^ dq 



H/VF 1 



Now, dynamic forward-scattering (with transfers q u)/vf) is not a perturba- 
tion anymore: on the contrary, the e dependence of ImS^ is dominated by for- 
ward scattering (the In |e|-term is larger than the "any-angle" e^-contribution 
). Correspondingly, the real part acquires a non-analytic term ReS oc e\e\, and 
the specific heat behaves as ^ 



The non-analytic T^-term in the specific heat has been observed in recent exper- 
iments on monolayers of He^ adsorbed on a soUd substrate [43]^. 

Finally, in ID the same power-counting argument leads to ImE^ cx |e| and 
ReE^ oc e In |£| '"^Correspondingly, the "correction" to the specific heat behaves 
as TlnT and is larger than the leading, T— term. This is the ultimate case of 
dynamic forward scattering, whose precursors we have already seen in higher 
dimensions 

** again, only processes with small total momentum contribute 

'if a term in C{T) does not fit your definition of non-analyticity, you have to recall that the 
right quantity to look at is the ratio C{T)/T. Analytic behavior corresponds to series C{T)/T = 
7 + 5T'^ + (tT* + . . . . whereas we have a T'^ In T and T terms as the leading order corrections to 
the Sommerfeld constant 7 for D = 3 and D = 2, correspondingly. 

'"Special care is required in ID as in the perturbation theory one gets a strong divergence in the 
self-energy corresponding to interactions of fermions of the same chirality (Fig. 8a,c). This point will 
be discussed in more detail in Section 2.3 (along with a weaker but nonetheless singularity in 2D). 
For now, let us focus on a regular part of the self-energy corresponding to the interaction of fermions 
of opposite chirality (Fig. 8b). 

'^Bosonization predicts that C(T) of a fermionic system is the same as that of ID bosons, which 
scales as T for ZP = 1 [10]. This is true only for spinless fermions, in which case bosonisation 
provides an asymptotically exact solution. For electrons with spins, the bosonized theory is of the 
sine-Gordon type with the non-Gaussian (cos<^) term coming from the backscattering of fermions of 
opposite spins. Even if this term is marginally irrelevant and flows down to zero at the lowest energies, 
at intermediate energies it results in a multiplicative InT factor in C(T) and a lnmax{g, T, H} 
correction to the spin susceptibility (where H is the magnetic field, and units are such that q, T, and 
H have the units of energy). The difference between the non-analyticities in D > 1 and D = 1 is 
that the former occurs already at the second order in the interaction, whereas the latter start only at 
third order. Naive power-counting breaks down in ID because the coefficient in front of T In T term 
in C(T) vanishes at the second order, and one has to go to third order. In the sine-Gordon model, the 
third order in the interaction is quite natural: indeed, one has to calculate the correlation function of 
the coscti term, which already contains two coupling constants; the third one occurs by expanding the 
exponent to leading (first) order. For more details, see [47], [48], [49]. 



C{T) = 72T - 02T\ 
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Even if the bare interaction is point-like, the effective one contains a long- 
range part at finite frequencies. Indeed, the non-analytic parts of S and C(T) 
come from the region of small q, and hence large distances. Already to the sec- 
ond order in U, the effective interaction tj = U^H^lo, q) is proportional to the 
dynamic polarization bubble of the electron gas, Il{uj,q). In all dimensions, 
Imll^ is universal and singular in q for \u!\ /vp ^ q ^ kp 

Imn^ (w, q) ~ z/n— 

VF\q\ 

Although the effective interaction is indeed screened at g ^ -and this is 
why the FL survives even if the bare interaction has a long-range tail-it has a 
slowly decaying tail in the intermediate range of q. In real space, U{r) behaves 
as uj/r^~^ at distances kp^ <C r <c vp/\ijj\. 

Thus, we have the same singular behavior of the bubble in all dimensions, and 
the results for the self-energy differ only because the phase volume q^ is more 
effective in suppressing the singularity in higher dimensions than in lower ones. 

There is one more special interval of q: q ~ 2kF , i-e., Kohn anomaly. Usu- 
ally, the Kohn anomaly is associated with the 2kF- non-analyticity of the static 
bubble, and its most familiar manifestation is the Friedel oscillation in electron 
density produced by a static impurity (discussed later on). Here, the static Kohn 
anomaly is of no interest for us as we are deahng with dynamic processes. How- 
ever, the dynamic bubble is also singular near 2kF- For example, in 2D, 

Imn^ {q » 2kF, to) oc ^ =9 {2kF - q) ■ 

^/kF {2kF - q) 

Because of the one-sided singularity in Imll^ near q = 2kF, the effective in- 
teraction oscillates and falls off as a power of r. By power counting, if a static 
Friedel oscillation falls off as sin 2kFr/r^, then the dynamic one behaves as 

ujsm2kFr 



Dynamic Kohn anomaly results in the same kind of non-analyticity in the self- 
energy (and thermodynamics) as the forward scattering. The "dangerous" range 
of q now is \q — 2fci?| ~ a;/fi;'-"dynanaic backscattering". It is remarkable that 
the non-analytic term in the self-energy is sensitive only to strictly forward or 
backscattering events, whereas processes with intermediate momentum transfers 
contribute only to analytic part of the self-energy. To see this, we perform the 
analysis of kinematics in the next section. 
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Fig. 6. Scattering processes responsible for divergent and/or non-analytic corrections to tlie self- 
energy in 2D. a) "Forward scattering"-an analog of the "g4"-process in ID. All four fermionic 
momenta are close to each other, b) Backscattering-an analog of the "(?2"-process in ID. The net 
momentum before and after colHsion is small. Initial momenta are close to final ones. Although the 
momentum transfer in such a process is small, we still refer to this process as "backscattering" (see 
the discussion in the main text). c)2kp— scattering. 



2.2. ID kinematics in higher dimensions 

The similarity between non-FL behavior in 1 D and non-analytic features in higher 
dimensions occurs already at the level of kinematics. Namely, one can make a 
rather strong statement: the non-analytic terms in the self-energy in higher di- 
mensions result from essentially ID scattering processes. Let's come back to 
self-energy diagram 3a. In general, integrations over fermionic momentum p and 
bosonic q are independent of each other: one can first integrate over {p, e), form- 
ing a bubble, and then integrate over (if, lo). Generically, p spans the entire Fermi 
surface. However, the non-analytic features in S come not from generic but very 
specific p which are close to either to k or to —k. 

Let's focus on the 2D case. The In |e| term results from the product of two 
-singularities: one is from the angular average of ImG and the other one 
from the dynamic, uj/vpq, part of the bubble. In Appendix Appendix A, it is 
shown that the Lojvpq singularity in the bubble comes from the region where p is 
almost perpendicular to q. Similarly, the angular averaging of ImG also pins the 
angle between k and q to almost 90°. 

ImG"^(e — w, fc — q) = — -Kb [e — lo — qvp coiiQ') ^ 

cose' = — ^i^e' 

vpq vpq 

As p and k are almost perpendicular to the same vector (q), they are either almost 
parallel or anti-parallel to each other. In terms of a symmetrized ("sunrise") self- 
energy (cf. Fig. 3), it means that either all three internal momenta are parallel to 

the external one or one of the internal one is parallel to the external whereas the 
other two are anti-parallel Thus we have three almost ID processes: 

'^In 3D, conditions p -L q and k J- q mean only that p and k he in the same plane. However, 
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• all four momenta (two initial and two final) are almost parallel to each other; 

• the total momentum of the fermionic pair is near zero, whereas the transferred 
momentum is small; 

• he total momentum of the fermionic pair is near zero, whereas the transferred 

momentum is near 2k p- 

These are precisely the same ID processes we are going to deal with in the 
next Section-the only difference is that in 2D, trajectories do have some angular 
spread, which is of order \uj\ /Ep. The first one is known as "(74" (meaning: all 
four momenta are in the same direction) and the other one as "(72" (meaning: two 
out of four momenta are in the same direction). Both of these processes are of 
the forward-scattering type as the transferred momentum is small. In ID, these 
processes correspond to scattering of fermions of same (54) or opposite chirality 
(52)- The last (2^^) process is known "g'l" in ID. 

It turns out that of these two processes, the 92- and 2kp- ones, are directly 
responsible for the In e behavior. The 54 -process leads to a mass-shell singu- 
larity in the self-energy both in ID and 2D, discussed in the next section, but does 
not affect the thermodynamics, so we will leave it for now. 

What about 2kF— scattering? Suppose electron k scatters into —k emitting 
an electron-hole pair of momentum 2k. In general, 2fc of the e-h pair may consist 
of any two fermionic momenta which differ by 2fc : p and p + 2k. But since 



2k 



2kF, the components of the e-h pair will be on the Fermi surface only if 



p « —k and p + 2k w k. Only in this case does the effective interaction (bubble) 
have a non-analytic form at finite frequency. Thus 2kF- scattering is also of the 
ID nature for £) > 1. 

What we have said above, can be summarized in the following pictorial way. 
Suppose we follow the trajectories of two fermions, as shown in Fig. 7. There 
are several types of scattering processes. First, there is "any-angle" scattering 
which, in our particular example, occurs at a third fermion whose trajectory is not 
shown. This scattering contributes regular, FL terms both to the self-energy and 
thermodynamics. Second, there are dynamic forward-scattering events, when 
q ^ \uj\ jvp. These are not ID processes, as fermionic trajectories enter the 
interaction region at an arbitrary angle to each other. In 3D, a third order in 
such processes results in the non-analytic behavior of C(T)-this is the origin 



it is still possible to show that for a closed diagram, e.g., thermodynamic potential, p and fc are 
either parallel or anti-parallel to each other. Hence, the non-analytic term in C (T) also comes from 
the ID processes. In addition, there are dynamic forward scattering events (marked with a star in 
Fig. 7) which, although not being ID in nature, do lead to a non- analydcity in 3D. Thus, the T'^ in T 
anomaly in C'(T) comes from both ID and non- ID processes [50] . The difference is that the former 
start already at the second order in the interaction whereas the latter occur only at the third order. In 
2D, the entire T— term in C{T) comes from the ID processes. 
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Fig. 7. Typical trajectories of two interacting fermions. Explosion: "any-angle" scattering at a third 
fermion (not sliown) wliicli leads to a regular (FL) contribution. Five-corner star: dynamic forward 
scattering q \lo\/vp. This process contributes to non-analyticity in 3D (to third order in the inter- 
action) but not in 2D. Four-corner star: ID dynamic forward and backscattering events, contributing 
to non-analyticities both in 3D and 2D. 



of the "paramagnon" anomaly in C (T) . In 2D, dynamic forward scattering does 
not lead to non-analyticity. Finally, there are processes, marked by"(7i", "52 
and "54", when electrons conspire to align their initial momenta so that they are 
either parallel or antiparallel to each other. These processes determine the non- 
analytic parts of Yj and thermodynamics in 2D (and also, formally, for D < 2.) 
A crossover between D > 1 and D = 1 occurs when all other processes but gi , 
f/2 , and (74 are eliminated by a geometrical constraint. 

We see that for non-analytic terms in the self-energy (and thermodynamics), 
large-angle scattering does not matter. Everything is determined by essentially 
ID processes. As a result, if the bare interaction has some q dependence, only 
two Fourier components matter: U (0) and L/(2fci?). For example, in 2D 

ImS^(£) cx [U^ (0) + {2kF) -U{0)U{2kF)]s^lii\s\; 

ReS^(w) cx [U^{0) + U'^{2kF)-U{0)U{2kF)]s\e\; 

C{T)/T = 7* - a [C/^ (0) + (2fcF) - ?7(0)t/(2/cF)] T; 

Xs{Q,T) = x:(0) + 6{/2(2fcf^)max{i;j.Q,T}; 

where a and b are coefficients. These perturbative results can be generalized for 
the Fermi-liquid case, when the interaction is not necessarily weak. Then the 
leading, analytic parts of C(T) and Xs are determined by the angular harmonics 
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of the Landau interaction function 

F{pJ)^F,{9)i + Fa{0)a-a', 

where 9 is the angle between p and p'. In particular, 

7* = 7o(l + (cos^J',)); 
^ o l + (cosgF.) 

where 70 and are the corresponding quantities for the Fermi gas. Because of 
the angular averaging, the FL part is rather insensitive to the details of the in- 
teraction. As generically Fg and Fa are regular functions of 6, the whole Fermi 
surface contributes to the FL renormalizations. Vertices U{0) and U{2kp), oc- 
curring in the perturbative expressions, are replaced by scattering amplitudes at 
angle 6 = n 

A{p,0)=A, {6)1 + A, {6) a -a', 
Beyond the perturbation theory [37], 

C{T)/T = r--a[Ali7:) + 3Al{7:)]T; 
Xs{Q,T) = x*s{0) + bAl{TT)max{vFQ,T}. 

Non-analytic parts are not subject to angular averaging and are sensitive to a 
detailed behavior of As^a near = tt". 

2.3. Infrared catastrophe 
2.3.1. ID 

By now, it is well-known that the FL breaks down in ID and an attempt to ap- 
ply the perturbation theory to ID problem results in singularities. Let's see what 
precisely goes wrong in ID. 1 begin with considering the interaction of fermions 
of opposite chirality, as in diagram Fig. 8b. Physically, a right-moving fermion 
emits (and then re-absorbs) left-moving quanta of density excitations (same for 
left-moving fermion emitting/absorbing right-moving quanta). Now, instead of 
the order-of-magnitude estimate (2.7), which is good in all dimensions but only 
for power-counting purposes, I am going to use an exact expression for the bub- 
ble, Eq. (B. 4), formed by left-moving fermions. On the Fermi surface (fc = kp, 

''The renormalization of the scattering amplitudes by the Cooper channel of the interaction results 
in additional In T-dependences of Ac,s(7r) 
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c) 

—►-J — ^ — ^ — — >■ 
+ + + + + 



Fig. 8. Self-energy in ID. ± refer to rigiit (left)-moving fermions. 



we have 

-ImE^. (e) ~ U'^vi duj J dqlmG^ {e-uj,k~q) Imli^ 

^ U'^vi j duj j dq6 {e — uj + vpq) {t^/vp) S {lj + vpq) 
- 9'\e\, 

where g = U/vp is the dimensionless coupUng constant. The corresponding real 
part behaves as e In \e\ . What we got is bad, as ImS^ scales with e in the same 
way as the energy of a free excitation above the Fermi level and ReE^ increases 
faster than e (which means that the effective mass depends on e as In |e|), but not 
too bad because, as long as g ^ 1, the breakdown of the quasi-particle picture 
occurs only at exponentially small energy scales: e < Ep ex-p{—l/g^). Now, 
let's look at scattering of fermions of the same chirality. This time, I choose to 
be away from the mass shell. 



dcj / dqS {e — uj — vp {k " q)) uj6 {ijj — vpq) (2.8) 



ImGj 



imni 



e^5 (e - vpk) 



(2.9) 
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It is not difficult to see that the full (complex) self-energy is simply 

^++ ■n+ - (2-10) 

On the mass shell (e = vpk) we have a strong-delta-function-singularity. This 

anomaly was discovered by Bychkov, Gor'kov, and Dzyaloshinskii back in the 
60s [52], who called it the "infrared catastrophe". Indeed, it is similar to an 
infrared catastrophe in QED, where an electron can emit an infinite number of 
soft photons. Likewise, since we have linearized the spectrum, a ID fermion 
can emit an infinite number of soft bosons: quanta of charge- and spin-density 
excitations. The point is that in ID there is a perfect match between momentum 
and energy conservations for a process of emission (or absorption) of a boson 
with energy and momentum related by uj = vpq '■ 

k' = k-q 

e' = e — u) = vpk — vpq = vp {k — q) . 

On the mass-shell, the energy and momentum conservations are equivalent. Imag- 
ine that you want to find a probability of certain scattering process using a Fermi 
Golden rule. Then you have a product of two 6— functions: one reflecting the 
momentum and other energy conservation. But if the arguments of the delta- 
functions are the same, you have an essential singularity: a square of the delta- 
function. As a result, the corresponding probability diverges. 

A pole in the self-energy [Eq. (2. 10)] indicates the non-perturbativeand specif- 
ically ID effect: spin-charge separation. Indeed, substituting Eq. (2.10) we get 
two poles corresponding to excitations propagating with velocities vp (1 ± g) 
(recall that 51 <C 1). This pecuUar feature is confirmed by an exact solution (see 
Section 3): already the 34— interaction leads to a spin-charge separation (but not 
to anomalous scahng). What we did not get quite right is that the velocities of 
both-spin- and charge-modes-are modified by the interactions. In fact,the exact 
solution shows that the velocity of the spin-mode remains equal to vp, whereas 
the velocity of the charge mode is modified. 

Obviously, there is no spin-charge separation for spinless electrons. Indeed, in 
this case diagram Fig. 8a does not have an additional factor of two as compared to 
Fig. 8c (but is still of opposite sign), so that the forward- scattering parts of these 
two diagrams cancel each other. As a result, there is no infrared catastrophe for 
spinless fermions. 

2.3.2. 2D 

What we considered in the previous section sounds like an essentially ID effect. 
However, a similar effect exists also in 2D (more generally, for 1 < D < 2). 
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This emphasizes once again that the difference between D = 1 and -D > 1 is not 
as dramatic as it seems. 

In 2D, the self-energy also diverges on the mass shell, if one linearizes the 
electron's spectrum, albeit the divergence is weaker than in ID-to second order, 
it is logarithmic The origin of the divergence can be traced back to the form 
of the polarization bubble at small momentum transfer, Eq. (A. 1). Integrating 
over the angle in 2D, we get 

Imn«(a;,q) = -(^) ^=^^=9 {vpq - \lo\) . (2.11) 

Imn^(w,g) has a square-root singularity at the boundary of the particle-hole 
continuum, i.e., at w = vpQ- (This is a threshold singularity of the van Hove 
type: the band of soft electron-hole pairs is terminated at = vpq, but the 
spectral weight of the pairs is peaked at the band edge). On the other hand, 
expanding ek+q in G^{e + u, k-l-g) as ^k+q = Cfe + ^fQ cos 6 and integrating 
over 9, we obtain another square-root singularity 



/ 



deimG^ = -27r 



{vpqf - {s + u-^kf 



-1/2 

(2.12) 



On the mass shell (oj = ^fe), the arguments of the square roots in Eqs. (2.1 1) and 
(2.12) coincide, and the integral over q diverges logarithmically. The resulting 
contribution to ImS^ diverges on the mass shell (s = ^k) [53, 54, 55], [34], [37] 

ImSf Je,fc)=- — -— In- 



where A = s — and u = mU /2t:. The process responsible for the log- 
singularity is the "34" process in Fig. 6. On the other hand, 52 and gi processes 
give a contribution which is finite on the mass shell 

(The divergence at e = — is spurious and is removed by going beyond the log- 
accuracy [34], [37].) We see therefore that the familiar form of the self-energy in 
2D [e^ In jej, seeRef. [56]] is valid only on the Fermi surface (^fe = 0) . The loga- 
rithmic singularity in ImS^ on the mass shell is eliminated by retaining the finite 
curvature of single-particle spectrum (which amounts to keeping the q^ /2m term 
in This brings in a new scale jEp. The emerging singularity in (2.3.2) 

'*In 3D, there is no mass-shell singularity to any order of the perturbation theory. 
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is regularized at |e — jEp and the In |£| behavior is restored. How- 

ever, higher orders diverge as power-laws and finite curvature does not help to 
regularize them. This means that-in contrast to 3D-the perturbation theory must 
be re-summed even for an infinitesimally weak interaction. Once this is done, the 
singularities are removed. Re-summation also helps to understand the reason for 
the problems in the perturbation theory. In fact, what we were trying to do was 
to take into account a non-perturbative effect-an interaction with the zero-sound 
mode-via a perturbation theory. Once all orders are re-summed, the zero-sound 
mode splits off the continuum boundary-now it is a propagating mode with ve- 
locity c> Vf- This splitting is what regularizes the divergences. The resulting 
state is essentially a FL: the leading term in E behaves as In |e| . However, 
some non-perturbative features remain: for example, the spectral function ex- 
hibits a second peak away from the mass shell corresponding to the emission of 
the zero-sound waves by fermions. A two-peak structure of the spectral function 
is reminiscent of the spin-charge separation, although we do not really have a 
spin-charge separation here: in contrast to the ID case, the spin-density collec- 
tive mode Ues within the continuum and is damped by the particle-hole pairs. 



3. Dzyaloshinskii-Larkin solution of the Tomonaga-Luttinger model 

3.1. Hamiltonian, anomalous commutators, and conservation laws 

In the Tomonaga-Luttinger model |57|,[58J one considers a system of ID spin- 
1/2 fermions with a linearized dispersion. Only forward scattering of left- and 
right-moving fermions is taken into account ((?2 and g4— processes), whereas 
backscattering is neglected. This last assumption means that the interaction po- 
tential is of sufficiently long-range, so that U (2^^) U {0) . [We will come 
back to this condition later.] Coupling between fermions of the same chirality 
(54) is assumed to be different from coupling between fermions of different chi- 
rality (g'2). If the original Hamiltonian contains only density-density interaction, 
then g2 = 94. A difference between 52 and 34 leads to an unphysical (within 
this model) current-current interaction. We will keep g2 7^ g4, however, at the 
intermediate steps of the calculations as it helps to elucidate certain points. At 
the end, one can make 52 equal to 54 without any penalty. In addition, in some 
physical situations, 52 7^ 54 • In what follows I will follow the original paper 
by Dzyaloshinskii and Larkin (DL) [59] and a paper by Metzner and di Castro 

"For example, Coulomb interaction between the electrons at the edges of a finite-width Hall bar 
(in the Integer Quantum Hall Effect regime) has this feature: electrons of the same chirality are 
situated on the same edge, whereas electrons of different chirality are on opposite edges; hence the 
matrix elements for the g2— and 514- interactions are different. 
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[60], where the Ward identity used by Dzyaloshinskii and Larkin is derived in a 
detailed way. 

The Hamihonian of the model is written as 

H = Ho + Hint ; 

-f^int = H2 + H4, 

where 

k.cr 

is the Hamiltonian of free fermions (± denote right/left moving fermions and a 
is the spin projection) and 

with 

P±,a = X! "^±,<T + 5) "i.-^ (^) • 

k 

To avoid additional complications, I assume that the interaction is spin-independent. 
To simplify the notations and to emphasize the similiirity between this model and 
QED, I will set i;f to unity in this section. 

Introducing the chiral charge- and spin densities as 

P± = P±A+P±,l' 
P± = P±A - P±,l' 
and total charge density and current as 

f = Pi -pi, 
the interaction part of the Hamiltonian reduces to 

Hint = E ^ + 54) / (<?) P' i-q) + 1(94- 92) f (g) 3' i-q) ■ (3.1) 

As we have already said, for = 94, the interaction is of a pure density-density 
type. Notice also that the spin density and current drop out of the Hamiltonian- 
this is to be expected for a spin-invariant interaction. To make a link with QED, 
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let us introduce Minkowski current ji^ with = 0, 1 so that = (=io) and 
= j'^ (=-.?! ) • Then the interaction can be written as a 4-product of Minkowski 
currents in a Lorentz-invariant form 



where 



ffoo = ^(52+54); 
511 = ^ (54 - 52) ; 

501 = 510 = 0. (3.2) 
In what follows, we will need the following anomalous commutators 

[p±,a- (q) , -f^o] = ±Qp±,a- {q) ; 
[p±,<T (q) , H2] = ^-^qPT,'T (q) ; 

[P±,<T {q) , Hi] = ±^qp±^„ {q) . 

The derivation of these commutation relations can be found in a number of stan- 
dard sources [61, 10] and I will not present it here. Adding up the conomutators, 
we get 

[p±,a, H] = [p±^^, H0 + H2 + Hi] 

_i_ _i_ ff2 c _i_ 54 c 

Adding up equations for spin-up and -down fermions, we obtain 



[pl,H] =±qpl±^qp^^±^qpl. 

IT TV 

Finally, adding up the ± components yields 

idtp' = [p',H]=Vcqf, (3.3) 

where 

= IH 

TT 

(recall that vp = !)• Eq. (3.3) is a continuity equation reflecting charge conser- 
vation. As if we did not have enough new notations, here is another one 

Q^ = {iV,q) 
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Fig. 9. a) Three-leg correlator K. b) Vertex part A. 



and 



In these notations and after a Fourier transform, the continuity equation can be 
written as 

Q^j" = 0. 

The same relation for free particles reads 

O^j^ = 0. 

3.2. Reducible and irreducible vertices 

Now, construct a mixed (fermion-boson) correlator 

K'^^„ {k, q\t, ti,t[) = -{Tf {q, t) a±,^ {k, h) a^^^ [k + q, h)), (3.4) 

where /x = 0, 1 and 




Fig. 10. Relation between vertices A and F. 



is an analog of the three-leg vertex in QED , except that in QED the "boson" 
is the /X— the component of the photon field 

QED : K^" = -{TA'\m). 

A diagrammatic representation of is a three-particle (one boson and two 
fermions) diagram (cf. Fig. 9a). 

The diagrams with self-energy insertions to solid lines simply renormalize the 
Green's functions. Absorbing these renormalizations, we single out the vertex 
part, re-writing as 

Rf" = G^A''. (3.5) 

Notice that there are as many vertex parts as there are bosonic degrees of free- 
dom. In (3+1) QED, A° is a scalar vertex and A''"^'^'^ are the components of 
the vector vertex. Diagrams representing A^ are shown in Fig. 9b. These se- 
ries can be re-arranged further by separating the photon-irreducible vertex part, 
r^. A photon-irreducible part is obtained by separating the corrections to the 
bosonic line, i.e., taking into account polarization. Vertices A** and F'' are re- 
lated via a kind of Dyson equation, which is simpler than the usual Dyson in 
a sense that there is no A^ on the right-hand-side. Diagrammatically, this rela- 
tion is represented by Fig. 10 where a shaded bubble is an exact (renormalized) 
current-current correlation function 

A^^^ {q,t) = -^{r {q,t)f (-g,0)). 

Algebraically, equation in Fig. 10 says 

A^,. = r^,. + A'^''5,.Arl<,. (3.6) 

(We remind the reader that indices ±, a simply specify the fermionic flavor which 
is not mixed in our approximation of forward- scattering and spin-independent 
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forces, so all relations are applicable to each individual flavor). The coupling 
constants relate currents to densities. According to Eqs. (3.1) and (3.2), den- 
sities couple to densities and currents to currents with no cross terms. Opening 
the matrix product in Eq. (3.6), we obtain 

Af = rf + A''%or° + ^'''5iir'. (3.7) 

3.3. Ward identities 

A Ward identity for vertex A^ is obtained by applying idt to in Eq. (3.4) 
and using the continuity equation (3.4).^^ Performing this operations and Fourier 
transforming in time, we obtain 

Q^i^f {K, Q) = Gi (K) -GiiK + Q), 

where i denotes the branch 

i = ±,(7. 

Recalling Eq. (3.5), we see that the Ward identity becomes 

Q^Af {K, Q) - (K + Q)- G;' (K) , (3.8) 

which is identical to a corresponding identity in QED. For those who hke to see 
things not masked by fancy notations, here is Eq. (3.8) in an expUcit form 

ujA° {s, k; CO, q) - v,qh\ {s, k; uj, q) = G^' (e + uj,k + q) - GJ^ (e, k) . (3.9) 

Notice that Eqs. (3. 8, 3. 9) contain renormalized velocity Vc- In what follows, we 
will actually need a Ward identity not for A^ but for the photon-irreducible ver- 
tex F'^. This one is obtained by deriving the continuity equation for 4-current 
correlation function A^" . To this end, one applies idt to A^"^ and uses continuity 
equation (3.3), which yields 

Q^A^-^ = -q5^,i. (3.11) 

TT 

'*When differentiating, recall that the T— product can be represented by step-functions in time 
which, upon differentiating, yield delta-functions 

'^To get this result, recall the form of the anomalous density-density commutator 

[?■'' (?),r(-?)] = 6'"'-?L, 

where e°° = e^^ = 0, = -e^" = 1. Now open the T- product in A°'^ and apply idt 
idtA°- {q, t) = (ift) (0 (t) / t) 0) + e i-t) f i-q, 0) f {q, t)> 

= ^Sit)lf iq,0) ,ri-q,0)]+VaqA^-' =2^S^,i+VaqA^-'. (3.10) 

V TT 

In 4-notations, (3.10) is equivalent to (3.11). 
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Now, we form a scalar product between and Eq. (3.7), using continuity 
equation for A^^" (3.11). This brings us to 

where 

=0 =2q/-K 

= ujV'i- Vc qT^+—- 

= l + (g4-S2)/7r 

Finally, the Ward identity for photon-irreducible vertex is 

Q^T>' = GT^{K + Q)-G-^{K). (3.12) 

It is remarkable that the left-hand- side of Eq. (3.12) contains the bare Fermi 
velocity (= 1) instead of the renormaUzed one. This is true even if we allowed 
for spin-dependent interaction in the Hamiltonian. 

It seems that we have not achieved much, as the conservation law was simply 
cast into a different form. However, in our ID problem with a linearized spec- 
trum a further progress can be made because the current and density (for given 
chirality) are just the same quantity (up to an overall factor of the Fermi velocity): 

i ±,(7 — ±i ±,0- 

Therefore, we have a closed relation between just one vertex and Green's func- 
tions. Suppressing the 4-vector index /x, we get the Ward identity for the density 
vertex 

G+\ (K + Q)- G+i (K) 



This is the identity that we need to proceed further with the Dzyaloshinskii- 
Larkin solution of the Tomonaga-Luttinger problem. Notice that (3.13) contains 
fully interacting Green's functions. 
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Fig. 1 1 . Dyson equation for the effective interaction. Solid line: Green's function of a right-moving 
fermion. Dashed line: Green's function of a left-moving fermion. Single wavy Une: bare interaction 
of fermions of the same chirality; spiral line: same for the fermions of opposite chirality. Double 
wavy and spiral lines represent the renormalized interactions. 



3.4. Effective interaction 

Effective interaction is obtained by collecting polarization corrections to the bare 
one. Diagrammatically, this procedure is described by the Dyson equation, rep- 
resented in Fig.ll. The interaction and polarization bubble are matrices with 
components 

where we used an obvious symmetry V++ = V = V-+. The Dyson 

equation in the matrix form reads 

V = Vo + VotW, 

or, in components, 

V++ = 54 + giXl+V++ + g2^-V+-\ 

v+- = <?2 + 52n_y++ + <74n_i/+_. (3.14) 

The bubbles in these equations 'xe fully renormalized ones, i.e., they are built on 
exact Green's functions and contain a vertex (hatched corner): 



r r dhrip 
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Now we use the Ward identity for V\ (3.13) to get 



18 



n± {uj,q)=-2i 



1 



'e + LO,k + q)\. (3.15) 



Eq. (3.15) looks exactly the same as a free bubble [cf. Eq. (B. 1)] except that 
it contains exact rather than free Green's functions. Because we managed to 
transform the product of two Green's functions into a difference, frequency inte- 
gration in Eq. (3.15) can be performed term by term yielding exact momentum 
distribution functions n± {k) and n± {k + q) : 



At the first glance, it seems that we have not achieved much so far Indeed, we 
traded one unknown quantity (n±) for another (n±). Both of them include the 
interaction to all orders and without any further simplification we are stuck. In 
fact, we have already made an important simplification: when specifying the 
model, we assumed only forward scattering. This means that the interaction 
is sufficiently long-range in real space so that backscattering can be neglected. 
Equivalently, in the momentum space it means that our interaction operates only 
in a narrow window of width qo near the Fermi points, ikp- Thus the states far 
away from the Fermi points are not affected by the interaction. The momentum 
integral in (3.16) comes from regions far away from the Fermi surface where 
unknown functions n± can be approximated by free Fermi steps. This approx- 
imation is good as long as qo <C kp- The solution is going to be exact only in 
a sense that there will be no constraints on the amplitude of the interaction (pa- 
rameters 52 and 514) but not its range. '^Now we understand better why the title of 
the paper by Dzyaloshinskii and Larkin [59] is "Correlation functions for a one- 
dimensional Fermi system with long-range interaction (Tomonaga model)"^*^. 

'^I skipped over a subtlety related to the infinitesimal imaginary parts iO+in the denominator. 
Works the same way. If you are unhappy with this, imagine that we work with Matsubara frequencies. 
Then there are no iO+s whatsoever. 

"in higher dimensions, we have a familiar problem of the Coulomb potential. Because it's a 
power-law potential, one cannot separate it into "amplitude" and "range". There is in fact a sin- 
gle dimensionless parameter, rs, which must be small for the perturbation theory-Random Phase 
Approximation-to work. Once Ts <g 1, we have two things: the screened potential is simultaneously 
weak and long-ranged. The Tomonaga-Luttinger model unties these two things: the interaction is 
assimied to be long-ranged but not necessarily weak. 

^"What seemed to be just a matter of mathematical convenience in the 70s, turns out to be quite a 
realistic case these days. If a wire of width a is located at distance d to the metallic gate, the Coulomb 
potential between electrons in the wire is screened by their images in the gate. Typically, d ^ a. A 
simple exercise in electrostatics shows that in this case U{0) is larger than C/(2fe^) by large factor 
In (d/a) [62]. 




(3.16) 
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With this simpUfication, the momentum integration proceeds in the same way 
as for free fermions (see Appendix Appendix B) with the result that the fully 
interacting bubbles are the same as free ones 

n±{uj,q)=nl{uj,q) = ±- . (3.17) 

-K u) — q + zO+sgnw 

This is a truly remarkable result which is a cornerstone for the DL solution 

Because our bubbles were effectively "Uberated" from the interaction effects, 
system (3.14) is equivalent to what we would have obtained from the Random 
Phase Approximation (RPA). It turns out that RPA is asymptotically exact in ID 
in the limit qo/kp — > 0. Solving the 2 by 2 system, we obtain for the effective 
interaction 

, . , . , ^ 54 (w + (?) + (g| - gl) q/iT 
y++ ij^, q) = {oj-q) 



where 



For gi = g2= g, 



9 9 

LU^ — q'^ 



= g ,.2 .,2.2 , .n^ - (3-18) 



^2 _ ^2^2 ^ 

3.5. Dyson equation for the Green's function 

The Dyson equation for right-moving fermions reads 

E+ {P) = i[ (P - Q) V++ (Q) r° (P, Q) . 

Diagrammatically, this equation is shown in Fig. 12. For Unear dispersion, 

Si {e,p) = eTP- {e,p) 
Substituting this relation back into the Dyson equations, we obtain 

f f duo dq 

{s-p)G+{e,p) = l+i / ——^G+{s,p)G+{s-oj,p-q)V++{uj,q)T°_{e,p-,uj,q). 

J J [2n) 

^'in QED, this statement is known as Furry theorem (W. H. Furry, 1937) 

^^Notice that as long as §4 7^ g2, the left-right symmetry is broken, i.e., the potential is not 
symmetric with respect to q ^ —q. 




Fig. 12. Dyson equation for the self-energy. 



Using the Ward identity (3.13) , we get 

(e-p+I]o)G+ {e,p) = ^+ij / ^^*^+ ^^'P^ i^-'^^P-Q) 

oj — q ^ 
J J f27r) ' w-q 



(27r) 
where 

const =i 



dwdq V++ (a>, q) 

{2nf -q 



. A constant term can always be absorbed into E, which simply results in a shift 
of the chemical potential. We are free to choose this shift in such a way that 
const=0, so that the Dyson equation reduces to 

ie-p)G,ie.p) = l + ^| ( ^,G, [e - ..p - q)Y±±^. (3.19) 

Notice that Eq. (3.19) is an integral equation with a difference kernel, which 
can be reduced to a differential equation for G. Before we demonstrate how it is 
done, let's have a brief look at a case when there is no coupling between left- and 
right-moving fermions: 52 = 0. In this case, 

T. / A (w - 1) (w - q) 
y++ q) = TT r^KT^ 



where 

w = 1 + 54/77. 
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Eq. (3.19) takes the form 



{e-p) G+ {e,p) = l + i{w-l) 




dujdq G+ {e — Lv,p — q) 



An w — wq + iO+ 



This equation is satisfied by the following function 



G+ {e,p) = 



1 



— p + iO+-^e — wp + iO+ 



(3.20) 



This is an example of a non-Fermi-Uquid behavior: the pole of a free G splits into 
the product of two branch cuts, one peaked on the mass shell of free fermions 
(e = p) and another one at the renormaUzed mass shell (e = wp). As left- and 
right movers are totally decoupled in this problem, the same result would have 
been obtained for two separate subsystems of left- and right movers. For exam- 
ple, Eq. (3.20) predicts that an edge state of an integer quantum Hall system is 
not a Fermi liquid, if spins are not yet polarized by the magnetic field [63]. The 
same procedure for a spinless system would give us a pole-like G with a renor- 
malized Fermi velocity. The non-Fermi-hquid behavior described by Eq. (3.20) 
is rather subtle: it exists only if both e and p are finite. In the Umiting case of 
p = Q (tunneling DoS) we are back to a free-fermion behavior G [e, 0) = e~^. 
Also, the e— integral of Eq. (3.19) gives a step-hke distribution function in mo- 
mentum space. The spectral function, however, is characteristically non-FL-like: 
instead of delta-function peak we have a whole region \p\ < |e| < w \p\ in which 
ImG is finite. At the edges of this interval ImG has square-root singularities. 

3.6. Solution for the case 52 = 54 

Substituting the effective interaction (3.18) into Dyson equation (3.19), we obtain 



Notice that the constant g is replaced by a momentum-dependent interaction, 
g {q) . The reason is that without such a replacement the integral diverges at the 
upper hmit. Here, the assumption of a cut-off in the interaction becomes impor- 
tant again. Transforming back to real time and space 




where 



u=^l + 2g/n. 
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we obtain the Dyson equation in a differential form 
where 



— + —]G{x,t)=P{x,t)G (x, t) - i5 {x)S{t) , (3.21) 



n:c,t) = ^JJ c^.c^,e-(^--),(,) ^,_- + ;.^^ . (3.22) 

The integral for P diverges if g is constant. To ensure convergence, we will 
approximate g{q) = ge"'*'/'^''. An actual form of the cut-off function is not 
important as long as we are interested in such times and spatial intervals such 
that x,t ':$> q^^ . The integral over uj is solved by closing the contour around the 
poles of the denominator ui = ±u \q\ (1 + iO"*") . For f > 0, we need to choose 
the one with Imw < 0. Doing so, we obtain 



/ 



^ . . . = . sgng + " ^»(qx-«kt|) 
27r 2u 



Solving the remaining q— integral, we obtain for P{x,t) 



A-KU \x — ut + i/qQ x + ut + i/qo 

For <: < 0, one needs to change go —qo in the last formula. 
The delta-function term can be viewed as a boundary condition 

G{x,0+)-G{x,0-) = -iS{x). (3.23) 

Once the function P {x, t) is known, Eq. (3.21) is trivially solved in terms of new 
variables r = x — t, s = x + t. For example, for t > 



G+ (r, t>0) = Go {x, t) /> (r) exp 



i [ ds'P{r,s') 

J r 



(3.24) 



where function /> (r) is determined by the analytic properties of G as a function 
of e. Substituting result for P (x, t) into Eq. (3.24), we get 

^ / . ^ n\ 1 ^ ^ +^ -f /- ^\ f x-t + i/qo\°''^^^^ f x-t-i/qo 

G+ {x, i > 0) = —Go {x, t) /> {x-t)[ , , ., — — — 'ri— 

2n \x — ut + i/qoJ \x + ut — i/qo 

where 

a = — . (3.25) 

8u 
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Formula for t < is obtained by choosing another function /< and replacing 
90 — > — 9o- Functions />,< are determined from the analytic properties. First of 
all, recall that 

Go {x,t) = — 1 — . 

We see that although Go is not an analytic function of t for any it is analytic 
for Ret > in the right lower quadrant (Imt < 0) and for Ret < in the upper 
left quadrant (Imt > 0). The interaction cannot change analytic properties of a 
Green's function hence we should expect the same properties to hold for full G. 

23 

From the boundary condition (3.23), it follows that 

f>{x) = f<{x) 
and / (0) = 0. 

Analyzing different factors in the formula for G, we see that only the term 
(a; — t =F Vso)" does not satisfy the required analyticity property. This term is 
eliminated by choosing function / {x) as 

/(x) = (goV + l)-". 
Finally, the result for G takes the form 



G+ {x,t) 



1 1 / a; - i + i7 ^ 



21^ X — t + isgnfO+ \x — ut + i^ 
1 

X 

[lo {x — ut + i'y) {x + ut — 



where 7 = sgnt/go- It seems somewhat redundant to keep two different damp- 
ing terms (isgnfO"'" and 7) in the same equation. However, these terms contain 
different physical scales. Indeed, isgntO"'" enters a free Green's function and 0"^ 
there has to be understood as the Umit of the inverse system size. On the other 
hand, 7 contains a cut-off of the interaction. Obviously, I7I 1/_L ^ O"*" for 
a realistic situation. The difference between the two cutoffs becomes important 

^'indeed, this property follows immediately from the Lehmann representation for G 

Gix,t) = -j^ |M^oPe'''"'^'e"'<^--^«'*, fort > 0; 

= |M„o|^e-'P-^e'(^--^o)*, fort < 0, 

where M^q are the matrix elements between the ground state and state u with energy Ei, > Eg. The 
required property simply follows from the condition for convergence of the sirai. 
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for the momentum distribution function and tunneling DoS, discussed in the next 
Section. 



3. 7. Physical properties 
3.7.1. Momentum distribution 

Having an exact form of the Green's function, we can now calculate the momen- 
tum distribution of, e.g., right-moving fermions: 



n+ 



/CSO 
-co 

ipx ^ 



I 

i 

'2^ 



dx 



dxe 



X + tO+ [q^x^ + If 
1 



11 r 

= o ^S^P / 

2 TT Jo 



ipx 



dx 



V mb (a;) 

X 



sin \p\ X 



X [q^x'^ + if 



We are interested in the behavior at p ^ (which means \p\ ^ qo). The final 
result for n+ {p) depends on whether a is larger or smaller than 1/2 [59, 64]. 
• For a < 1/2 ("weak interaction"), one cannot expand sin pa; in x because the 
resulting integral diverges at a; = oo. Instead, rescale px ^ y 



b) = ^ - ^ 



dy 



smy 



and neglect 1 in the denominator. This gives 



{qo/p) 2/^ + 1 



sgnp 



(3.26) 



where 



smTra 



r(-2a). 



Notice that n+ (p) is finite (= 1/2) at p = 0, although its derivative is singular. 
We should be able to recover the Fermi-gas step at p = by setting a = in 
(3.26). Indeed, 

lim Ci = a = — 

a^o -2a 2 



and 



n{p) = 



1 — sgnp 
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which is just the Fermi-gas result. Notice also that there is nothing special about 
the limit a — > ^. Indeed, constant C\ has a regular expansion in a 

1 

Ci = -- -70; + ..., 

where7 = 0.577. . . andfactor /qof" can be expanded for finite p and small 
a as 

(bl /go)'" = l + 2aln|p| /go. 
To leading order in a, we obtain 

{P) = \- sgnpi [1 + 2a In \p\ /qo] = no (p) - asgnpln \p\ /qo, 

which is a perfectly regular in a (but logarithmically divergent at p ^ 0) be- 
havior. Once again, it is not surprising: despite the fact that the results for a ID 
system differ dramatically from that for the Fermi gas, they are still perturbative, 
i.e., analytic, in the coupUng constant. 

• For a > 1/2 ("strong interaction"), it is safe to expand sin pa; and the result 

is ^ 

(P) = 2 ~ C'ap/g'o, 

where 

1 r(a-l/2) 
' 2^ T{a) • 

In this case, no remains of a jump at the Fermi point is present in n+ (p) which 
is a regular, linear function near p = 0. 

• Finally, a = 1/2 is a special case, where expansion in p results in a log- 
divergent integral. To log-accuracy 

n+ (p) = In — . 

In general, n (p) is some hypergeometric function ofp/qo which decays rapidly 
forp ^ qo and approaches 1 forp <C —qo- A posteriori, this justifies the replace- 
ment of exact n (p) by its free form in the Dyson equation. 

3.7.2. Tunneling density of states 
Now we turn to the tunneling DoS 

N{e) = -ilmG^ (e,a; = 0). 

TT 

^^contrary to some statements in the literature. 
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Recalling that [23] 



we see that 



G(e),fore > 0; 
G* (e) , for e < 0, 



ImG-" (e, 0) = sgnelmG (e, 0) 



and 



N {e) = —— SgnelmG (0,£) = ——sgne 

TT TT 



1 

sgne 

TT 



j dte'^*G (0, t)- j dte-'^*G* (0, t) 
j dte''* {G (0, t) - G* (0, -t)} 



For f — > oo, 



and 



G(0,i) 



const 



l+2a 



(-0 

G{0,t)-G* i-t) 



is an odd function of t. Thus 



A^'fe) = — — sgnelmG (0,£) = — — sgne— 

TT TT 2i 



= sgne 

TT 



sin £t 







y dte'^*G (0, t)- j dte-'^^G* (0, i) 
/ dtsin£t{G(0,f) - G* (0,-t)} oc sg 
The integral is obviously convergent for a < 1/2. In this case, 

Ns (£) oc |£|'" , 

which means that the local tunneling DoS is suppressed at the Fermi level. Actu- 
ally, the exponent for a > 1/2 is the same, however, the prefactor is a different 
function of a [65]. 

The DoS in Eq. (3.7.2) with exponent 2a, where a is given by Eq. (3.25) 
corresponds to tunneling into the "bulk" of a ID system, i.e., when the tunnehng 
contact (with a tip of an STM or another carbon nanotube crossing the first one) 
is far away from its ends. In the next Section, we will analyze tunneling into an 
edge of a ID conductor, which is characterized by a different exponent, a' . 
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4. Renormalization group for interacting fermions 

The Tomonaga-Luttinger model can be solved exactly as it was done in the previ- 
ous Section- only in the absence of backscattering. Backscattering can be treated 
via the Renormalization Group (RG) procedure. This treatment is standard by 
now and discussed in a number of sources [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. For the 
sake of completeness, I present here a short derivation of the RG equations. A 
reader familiar with the procedure can skip this Section and go directly to Sec. 5, 
where these equations will be used in the context of a single-impurity problem. 

An exact solution of the previous Section is parameterized by two coupling 
constants, g2 and g^, which are equal to their bare values. In the RG language, 
it means that these couplings do not flow. Let's see if this is indeed the case. In 
what follows, I will neglect the 94— processes, as their effect on the flow of other 
couplings is trivial, and, for the sake of simplicity, consider a spin-independent 
interaction. To second order, the renormalization of the 52— couphng is ac- 
counted for by two diagrams: diagrams a) and b) of Fig. 13. 

Diagram a) is a correction to 52 in the particle-particle channel. The correc- 
tion to g2 is given by 



52 ) = — — 2 j ^1 j d<^G+ + ^w, fci -I- q) G- {ie2 - ico, k2- q) . 



Without a loss of generality, one can choose all momenta to be on the Fermi 
"surface": ki = k2 = = = 0. Choose q > (the other choice q < will 
simply double the result) 





Adding the result up with the (identical) g < contribution, we find 
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Diagram b) is a correction to 32 in the particle-hole channel: 

^ (27r) J J 

^ 5l 1 1 



dq / dw- 



(27r)^ io J i (£1 + w) - g i (£4 + w) + g 

T292 I Q- — :r—v^ = -—In 



(27r) Jo £1 - £4 + w + 2iq Aw e\ - £4 
As in the previous case, the final result is: 



■In- 
27r £1 — £4 



If we sum only the Cooper ladders, adding up more vertical interaction Unes 
to diagram a), the fuU vertex becomes 



1 + 52 In 



ei+£2 



(To keep track of the signs, one needs to recall that in Matsubara frequencies each 
interaction line comes with the minus sign from the expansion of the S— matrix). 
The resulting vertex blows up for attractive interaction {gi < 0) as £1 + £2 — > 0, 
which is nothing more than a Cooper instability. 

Likewise, untwisting the crossed lines in diagram b) and adding more interac- 
tion Unes, we get the particle-hole vertex 

r - g2 

l-52ln^" 



This vertex has an instability for repulsive interaction (52 > 0). In fact, none of 
these instabiUties occur. To see this, add up the results of diagrams a) and b) 



)a+6 27r 



512 , iA , iA 



In In - 

£1 -|- £2 £1 — £4 



27r £1 -I- £2 ' 



In the RG, one changes the cut-off and follow the corresponding evolution of the 
couplings. As the cut-off dependence cancelled out in the result for (g^^ ) , 

V / a+b 

coupling 92 remains invariant under the RG flow. 

Backscattering generates additional diagrams: diagrams c)-f) in Fig. 13. 
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Fig. 13. Second order diagrams for couplings g2 (solid wavy line) and g\ (dashed wavy line). Straight 
solid and dashed lines correspond to Green's functions of right- and left moving fermions, corre- 
spondingly. 
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Diagram c) describes repeated backscattering in the particle-particle channel, 
which is equivalent to forward scattering. Therefore, this diagram gives a correc- 
tion to coupling. Using the relation between G±, i.e., G± = — (Gip)* , we 
find 



(92^) = f dq f (hjG-{isi+iiO,q)G+{iSi + iuj,q) 

^ 'c (27r) J J 

J dq J du)G+ [iei + iuj, q) G- {154 + iuj, q) 



{2n) 
{2nf 



The last integral is the same as for (q^^^ ■ Thus, 



2-77 £1 -h £2 



The rest of the diagrams provide corrections to g\ . 

Diagram dl) is the same as diagram a) except for the prefactor being equal 
to gxgi : 



(2; 
91 



.9152 



'dl 27r e\ + £2 

Diagram d2) is a complex-conjugate of diagram dl). The sum of diagrams 
dl) and d2) is equal to 



(^^ ) dV-Vdl 



9192 , iA 5152 , -iA 
■ m : h — — m ■ 



27r £1 -I- £2 2n e\ -\- £2 



TT £1+62 

Diagrams e) is a polarization correction to the bare coupling: 

(51^^)^= fi(?n2fep (W = £1 - £2,g = 0). 

fermionic loop 

Using Eq. (B. 8), we obtain 

iV., A_ 

ki - £21 



/ (2)\ 2i ^ 



where A^^ is the degeneracy factor (=2 for spin 1/2 fermions, occupying a single 
valley in the momentum space). 
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Diagram fl) is the same as the bubble insertion, except for no minus sign, no 
degeneracy factor {N^) factor, and the overall coefficient is gig2: 



V 7/1 277 \e\-. 



A 

'/I 277 |£l-£2| 
Diagram f2) is equal to f 1). Their sum 

1 . A 



V //1 + /2 TT 



|£l - £21 

Collecting all contributions together, we obtain 

, «A gl , lA 5? -iA 
r2 = 52-;^ In TT" In ■ 



27r £i + £2 27r £1 — £4 27r £1 + £2 ' 



cancel out in the RG s 



sense 



5ifl2,„ A TV, 2,„ A , 1„ „ A 



Ti = 51 In ■ —g-^ In + -5152 In 



£l+£2 27r |£l-£2| TT |£l-£2| 

Second and fourth terms in Ti also cancel out in the RG sense. Changing the 
cut-off from A to A + A, we obtain two differential equations 



dl 27r' 



(4.1) 



where Z = In A. We see that a quantity 

f = r2 - -^Fi = const = 32 - -^51- (4.3) 

is invariant under RG flow, therefore its value can be obtained by substituting the 
bare values of the couphng constants (32 and gi) into (4.3). The RG-invariant 
combination is then 

f = 52-^5i- (4.4) 
For spinless electrons {Ng = 1), 

f = T2-Vi = U{Q)-U{2kF). 
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This last result can be understood just in terms of the Pauli principle. Indeed, 
the anti-symmetrized vertex for spinless electrons is obtained by switching the 
outgoing legs of the diagram (pi,p2 P^iPa)- To first order, 

r(pi,P2;P3,P4) = U {pi -P3)-U (pi-pi) . 

Choosing P3 = pi — q andp4 = P2 + Q, we obtain [recall that U {q) = U {—q)] 

r(pi,P2|g) = U (q) -U{pi-p2- q). 

One of the incoming fermions is a right mover (pi = pf) and the other one is a 
left mover (p2 = —Pf)- As g is small compared to pf, we obtain 

T{pi,p2\q) = U{{))-U{2kF). 

In fact, for spinless electrons 52 and g\ processes are indistinguishable^^ as we 
do not know whether the right-moving electron in the final state is a right-mover 
of the initial state, which experienced forward scattering, or the left-mover of the 
initial state, which experienced backscattering. A proper way to treat the case of 
spinless fermions is to include backscattering into Dzyaloshinskii-Larkin scheme 
from the very beginning, re- write the Hamiltonian in terms of forward scattering 
with invariant couphng f , and proceed with the solution. All the results wiU then 
be expressed in terms of F rather than of 52- 
Solving the equation for Fi, gives on scale e 

^ (5ir^+WnA/e- ^"-'^ 

At low energies, Fi renormalizes to zero (T\ = Ti{l — 00) = 0), if the inter- 
action is repulsive, and blows up at e = Aexp (— l/|(7i|), if the interaction is 
attractive. Coupling r2 also flows to a new value which can be read off from 
Eq. (4.4) 

1 



^2=92- —91- 



Roughly speaking, gi is not important for repulsive interaction as the effective 
low-energy theory will look like a theory with forward scattering only. This does 

^^That does not mean that backscattering is unimportant! It comes with a different scattering 
amplitude U (2k p) . In fact, it is only backscattering which guarantees that the Pauli principle is 
satisfied, namely, for a contact interaction, when U (0) = U {2k p) , we must get back to a Fermi 
gas as fermions are not allowed to occupy the same position in space and hence they cannot interact 
via contact forced. Our invariant combination U (0) — U {2k p ) obviously satisfies this criterion. We 
will see that bosonization does have a problem with respecting the Pauli principle, and it takes some 
effort to recover it. 
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not really mean, however, that one can consider a fixed point as a new problem 
in which backscattering is absent, and apply our exact solution to this problem. 
Instead, one should calculate observables, derive the RG equations for flows, and 
use current values of coupling constants in these RG equations. An example of 
this procedure will be given in the next Section, where we will see that the flow 
of Fi provides additional renormalization of the transmission coefficient in an 
interacting system. 

Assigning different coupling constants to the interaction of fermions of par- 
allel (gi||) and anti-parallel {gi±) spins, one could see that the coupling which 
diverges for attractive interaction is in fact gi±. This clarifies the nature of the 
gap that RG hints at (in fact, a perturbative RG can at most just give a hint): it is a 
spin gap. This becomes obvious in the bosonization technique, as the instability 
occurs in the spin-sector of the theory. An exact solution by Luther and Emery 
[66] for a special case of attractive interaction confirms this prediction. 

5. Single impurity in a ID system: scattering theory for interacting fermions 

A single impurity or tunneling barrier placed in a ID Fermi gas reduces the con- 
ductance from its universal value-e-^/ h per spin orientation-to 

a = iV.^|io|', (5.1) 

where to is the transmission amplitude. The interaction renormalizes the bare 
transmission amplitude. As a result, the conductance depends on the charac- 
teristic energy scale (temperature or appUed bias), which is observed as a zero- 
bias anomaly in tunneUng. This effect is not really a unique property of ID 
: in higher dimensions, zero-bias anomalies in both dirty and clean (ballistic) 
regimes [67, 12, 13, 14] as well as the interaction correction to the conductiv- 
ity [67, 15], stem from the same physics, namely, scattering of electrons from 
Friedel oscillations produced by tunneling barriers or impurities. ID is special 
in the magnitude of the effect: the conductance varies significantly already on 
the energy scale comparable to the Fermi energy, whereas in higher dimensions 
the effect of the interaction is either small at all energies or becomes significant 
only at low energies (below some scale which is much smaller than Ep &i long 
as the parameter kpl, where I is the elastic mean free path, is large. The ID 
zero-bias anomaly is described quite simply in a bosonized language [68], which 
does not require the interaction to be weak. We will use this description in Sec. 6. 
However, in this Section I will choose another description-via the scattering the- 
ory for fermions rather than bosons-developed by Matveev, Yue, and Glazman 
[11]. Although this approach is perturbative in the interaction, it elucidates the 
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underlying mechanism of the zero-bias anomaly and allows for an extension to 
higher-dimensional case (which was done for the case of tunneling in Ref.[12] 
and transport in Ref.[15]). 



5.7. First-order interaction correction to the transmission coefficient 

In this section we consider a ID system of spinless fermions with a tunneling 
barrier located at a; = [1 1]. For the sake of simplicity, I assume that the barrier 
is symmetric, so that transmission and reflection amphtude for the waves coming 
from the left and right are the same. Also, I assume that e-e interaction is present 
only to the right of the barrier, whereas to the left we have a Fermi gas. Such 
a situation models a setup when a tunneUng contact separates a ID interacting 
system (quantum wire or carbon nanotube) and a "good metal", where one can 
be neglect the interaction. We also assume that the interaction potential U {x) 
is sufficiently short-ranged, so that U (0) is finite and one can neglect over-the- 
barrier interaction. However, U (0) ^ U {^kp) (otherwise, spinless electrons do 
not interact at aU^^). 

The wave function of the free problem for a right-moving state is: 

^l{x) = -l={e'>'-+roe-"^''),x<0; (5.2) 

= -^toe'^\x>0. 
VL 

For a left-moving state: 



toe-''''',x <0. (5.3) 



Here k = \/2mE > 0. To begin with, we consider a high barrier: |to| 1, ''o ~ 
— 1 . Then the free wavefunction reduces to 

V'fc (x) = — ^ sin fca;, a; < (incoming from the left+reflected); (5.4) 

= ^ioe''''',a; > (transmitted left ^ right); (5.5) 



*For a contact potential [which leads to f7 (0) = C/ (2feji')], the four-fermion interaction for the 

iless case reduces to 
the interaction is absent. 



spinless case reduces to [*t (0)]^ *2 (q) . By Pauli principle, [*t (o)]^ = 'J'^ (q) = q, so that 




Fig. 14. Correction to the Green's function: exact with respect to the barrier and tirst order in the 
interaction. 



7p\{x) = —j=tQe *''^,a;<0 (transmitted right -> left); (5.6) 

2i 

= -= sin kx, X > (incoming from the right+reflected)(5.7) 



The barrier causes the Friedel oscillation in the electron density on both sides of 
the barrier. The interaction is treated perturbatively, via finding the corrections to 
the transmission coefficient due to additional scattering at the potential produced 
by the Friedel oscillation. Diagrammatically, the corrections to the Green's func- 
tion are described by the diagrams in Fig. 14, where a) represents the Hartree and 
b) the exchange (Fock) contributions, correspondingly. Compared to the text- 
book case, though, the solid lines in these diagrams are the Green's functions 
composed of the exact eigenstates in the presence of the barrier (but no interac- 
tion). Because the barrier breaks translational invariance, these Green's functions 
are not translationally invariant as well. I emphasized this fact by drawing the di- 
agrams in real space, as opposed to the momentum -space representation. Notice 
also that the Hartree diagram is usually discarded in textbooks because the bub- 
ble there corresponds to the total charge density (density of electrons minus that 
of ions), which is equal to zero in a translationally invariant and neutral system. 
However, what we have in our case is the local density of electrons at some dis- 
tance from the barrier Friedel oscillation is a relatively short-range phenomenon 
(the period of the oscillation is comparable to the electron wavelength), and it is 
possible to violate the charge neutrality locally on such a scale. As a result, the 
Hartree correction is not zero. 

To first-order in the interaction, an equivalent way of solving the problem is 
to find a correction to the wave-function, rather than the Green's function, in 
the Hartree-Fock method. The electron wave-function which includes both the 
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barrier potential and the electron-electron interaction is 

Mx) = i^k{x)+ / dx'G>{x,x',E) 
Jo 

/■oo 

X / dx"[V„{x")S{x' -x") + V,,{x',x")]iljt{x"), (5.8) 
Jo 

where Gq is the Green's function of free electrons on the right semi-line, E 
is the full energy of an electron, Vh and V^x are the Hartree and the exchange 
potentials. The Hartree potential is 

^„(., = /^'a(»-.',M.'), (5.9) 

where Sn{x) = n{x) ~ uq is the deviation of the electron density from its uni- 
form value (in the absence of the potential) and U (x) is the interaction potential. 
Hartree interaction is a direct interaction with the modulation of the electron den- 
sity by the Friedel oscillation. For a high barrier, which is essentially equivalent 
to a hard-wall boundary condition, the electron density is 

nix) = Aj^ ^-^(M=-o(l-^^j- (5.10) 
8n{x) = -sm{2kFx) /2'KX, (5.11) 
where no = Aif/tt is the density of electrons. Then, 

/ N 1 1 iTT, ,.sm2kFx' ,^ -.^^ 

Vh{x) = -— dx'U{x-x') . (5.12) 

Notice that although the bare interaction is short-range, the effective interaction 
has a slowly-decaying tail due to the Friedel oscillation. (The integral goes over 
only for positive values of x' because electrons interact only there.) 
The exchange potential is equal to 



Vex{x,x) = -U{x-x')[j^ ^[^^(^')] ^kix) 



(5.13) 



Since we assumed that electrons interact only if they are located to the right of 
the barrier, the integral in (5.8) runs only over x, x' > and the Green's function 
is a Green's function on a semi-line. The wave-function in (5.8) needs to be 
evaluated at a; oo, which means that we will only need an asymptotic form of 
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the Green's function far away from the barrier. This form is constructed by the 
method of images 

C^ix^x'^E) = Gq{x,x',E) - Go{x,-x',E), (5.14) 

where 

Go {x,x',E) = J-e'H^-^'l 

IVk 

is the free Green's function on a Une with k = \phnE and Vk = k/m. Coordi- 
nate x' is confined to the barrier, whereas a; ^ oo, thus x> x' and 

G'^{x,x',E) = - —sin{kx')e'''\ 
5.1.1. Hartree interaction 

Our goal is to present the correction to the wave-function for electrons going 
from a; < to a; > in the form 

Vfe-Vfe = ^^ie^'", (5.15) 

where 5t is the interaction correction to the transmission coefficient. Substituting 
(5.15) into (5.8), we obtain for the Hartree contribution to t 



St^ 2 f°° 

— = / da; sin kxe'^'^'^VH {x) . 

to vp Jo 



to Jo 

where one can replace Vk vp iJ^ all non-oscillatory factors. For a delta- 
function potential, U {x) = U6 (x) 

U sm2kFX 

Vh{x)=-- . (5.16) 

2lT X 

However, the S— function potential is not good enough for us, because the Hartree 
and exchange contributions cancel each other for this case. Friedel oscillation 
arises due to backscattering. With a httle more effort, one can show that U in the 
last formula is replaced by U {2kF) 

, , U {2kF) sm2kFX 
^= 2^ 



^^Notice that the sign of the Hartree interaction is attractive near the barrier (assuming the sign of 
the e-e interaction is repulsive at 2fcf): for x ^ 0, Vh {x) —U {2k p) kp/ir. The reason is that 
the depletion of electron density near the barrier means that the positive background is uncompen- 
sated. As a result, electrons are attracted to the barrier and transmission is enhanced by the Hartree 
interaction. 
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Substituting this into yields 



where 



to ttvf Jo X 

U{2kF) f°° J i / 2ikx 1 , ■ , ^ sm2kFX 
> -J- I g^'i-i _ I regular correction to Imt I 



nvp 



f dx^ (e 
Jo 2^ V 



U{2kF) /■°"^^J_^2ifexSin2fcFa; _ [/(2A;f) 



ttvf Jo 2i X 2'kvf 

f°° , /' . ^, -1 r,, . , • \ sm2kFX 

: / ax I sm 2kx + i cos2fca; yet another regular correction 

U(2kF) f°° , . „, sin2A;Fa; 
' aa;sin2Ka;- 



2ttvf Jo ^ 

U{2kF) k + kF , kp 

1^ \k^\ \k-kF[ 



"2fc, 



Attvf 



and gi = U (2A;f) ■ In deriving the final result, all terms regular in the Umit 
k ^ kp were discarded. 



5.7.2. Exchange 

Now both X and x' > 0. We need to select the largest wave-function, i.e., such 
that does not involve a small transmitted component. Obviously, this is only 
possible for A: < (second term in (5.13)) and j, , given by (5.7). Substituting 
the free wave-functions into the equation for the exchange interaction, we get 

Vex{x, x') = -U{x - x')p{x, x'), (5.17) 

where the ID density-matrix is 

p{x,x') = 4/ — sin(fca;) sin(fca;') (5.18) 

Jo 277 

= 2 —[cosk{x-x')-cosk{x + x')] (5.19) 

Jo 27r 



sin/cF(a; -|- x') 
7r(x + x') 



(5.20) 
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where . . . stand for the term which depends on a; — a;' . This term does not lead to 
the log-divergence mSt and will be dropped . For a; = a;' , we get the correction 
to the density 5n (a;) , as we should. 

Correction to the transmission coefficient 

5t,,/to = — r dx' r dx"U {x' - x") sin kx'e^'^^" \ """^ . 

t^vf Jo Jo x' + x" 

(5.21) 

After a little manipulation with trigonometric functions, which involves dropping 
of the terms depending only on a; — a;', we arrive at 

- ' 'Sui,) r '-E± (5.22) 

X. 



to Att^vf Jo q 







where 



x{sin2(fc — /cf + q)x+ — sin2(A; — Uf — <l)x+}, (5.23) 



x+ = (5.24) 



Integral over x+ provides a lower cut-off for the q— integral 



— -{sm2{k-kF + q)x+-sm2{k-kF-q)x+} (5.25) 
^+ 

TT TT 

= -sgn(g + k-kF) + -^sgn{q -k + kF) (5.26) 

= 7761(9-1^-^^1). (5.27) 

Now 

^ = - 1 p ,5.28, 

to ^-KVF J\k-kF\ Q 

As U (q) is regular at g — > one can take U {q) out of the integral at g = 
(denoting U (0) = (72) 

Stex 1 r° dq qo 

-92 — = -tto In- 



to 'iTTVF J\k-kF\ 1 \k-kF\ 



^^Notice that the important part of the exchange potential is repulsive near the barrier This means 
that electrons are repelled from the barrier and transmission is suppressed. 

^'if U (q) has a strong dependence on q for q ^ (which is the case for a bare Coulomb potential 
U (q) oc in q), this dependence affects the resulting dependence of the transmission coefficient on 
energy \k — kp \ , i.e., on the temperature and/or bias. Instead of a familiar power-law scahng of the 
tunneling conductance for the short-range interaction, the conductance falls off with energy faster 
than any power law for the bare Coulomb potential. 
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Combining the exchange and Hartree corrections together (in doing so, we choose 
the smallest upper cut-off for the log which we assume to be the inverse interac- 
tion range, go)) we get 

Jt = -Wln |, | , (5.29) 
\k — kp] 

where 

a' = a'n — OL-ih^ = — — asymmetric geometry. (5.30) 

It can be shown in a similar manner that if we had interacting regions on both 
sides of the barrier, the result for ol would be double of that in Eq. (5.30). 

a' = ag — a2fep = ^ — symmetric geometry. (5.31) 

The sign of the correction to t depends on the sign of 92 ^ .9i = ?7 (0) — [/ {2kp) . 
Notice that transmission is enhanced, if U (2fci?) > J7(0). Usually, this behavior 
is associated with attraction. We see, however, that even if the interaction is 
repulsive at all q but U (^kp) > U{0), it works effectively as an attraction. 
The case U (2/0^) > U (0) is not a very realistic one, at least not in a situation 
when electrons interact only among themselves. Other degrees of freedom, e.g., 
phonons, must be involved to give a preference to 2kF— scattering. 

5.2. Renormalization group 

It is tempting to think that the first-order in interaction correction to to in Eq. (5 .30) 
is just an expansion of the scaling form tcx |fc — .A poor-man RG indeed 
shows that this is the case. Near the Fermi level, k — kp = {E — Ep) jvp = 
e/vp so that the first-order correction to t is 

ti =to [i- — a m -J— - 



where Wq — qovp is the effective bandwidth. The meaning of this bandwidth 
is that the states at ±Wo from the Fermi level (=0) are not affected by the inter- 
action. For |£| = Wo, ti = to- Suppose that we want to reduce to bandwidth 
Wo^Wi< Wo and find t at |e| = Wi 

t^=to \^l-a\n — 

It is of crucial importance here that coefficient a' (which will become the tun- 
neling exponent in the scaling form we are about to get) is proportional to the 
RG-invariant combination [/ (0) — [/ {2kp) =52 — 51 for spinless electrons. 



Fundamental aspects of electron correlations and quantum transport 



59 



This means that a' is to be treated as a constant under the RG flow. Repeating 
this procedure using t, found at the previous stage instead of a bare to, n times, 
we get 

tn+\ = 1 - a In 



Wn+l. 

The renormalization process is to be stopped when the bandwidth coincides with 
the physical energy |£| , at which t is measured. In the continuum limit (f„+i — 
tn = dt; Wn+i = Wn — dW), this equation reduces to a differential one 

dt _ ,dW 

Integrating from t (e) to to (and, correspondingly, from = |£| to = Wo), 
we obtain 

t{e)=to{\e\/Wor' . 



5.3. Electrons with spins 

Now let's introduce the spin. The effect will be more interesting than just mul- 
tiplying the result for the tunneling conductance by a factor of two (which is 
all what happens for non-interacting electrons.) To keep things general, I wiU 
assume an arbitrary "spin" (which may involve other degrees of freedom) degen- 
eracy Ns and put A^^ = 2 at the end. In this section we will exploit the result of 
Sec.4 stating the backscattering amplitude flows under RG. This flow affects the 
renormalization of the transmission coefficient at low energies. 

Repeating the steps for the first-order correction to t for the case of electrons 
with spin is straightforward: one just has to recall that the Hartree correction is 
multiplied by Ns (as the polarization bubble involves summation over all isospin 
components, it is simply multiplied by a factor of Ns). On the contrary, the ex- 
change interaction is possible only between electrons of the same spin, so there 
are Ns identical exchange potentials for every spin component. I am going to 
discuss the strong barrier case first in the symmetric geometry. Then, taking into 
account what we have just said about the factor of Ns, we can replace the result 
for spinless electrons (5.30) by 

a ^ a =ao- Nsa2kp = — • (5.32) 

(and similarly for the symmetric geometry of the tunneling experiment). Corre- 
spondingly, the correction to the transmission coefficient (for a given spin pro- 
jection) changes to 

t„ = to{l-a'lnW/\e\). 
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The tunneling conductance is found from the Landauer formula 



,2 



where, as the barrier is spin-invariant, the sum simply amounts to multiplying the 
result for a given spin component by Ns ■ Now, the result in Eq. (5 .32) seems to be 
interesting, as the 2fci? contribution gets a boost. If NgU (2A;f) > U (0), we have 
in increase of the barrier transparency. It does not seem too hard to satisfy this 
condition. For example, it is satisfied already for the delta-function potential 
and Ns = 2. However, as opposed to the spinless case, a' is not an RG-invariant 
but flows under renormalizations. Let's split a' into an RG-invariant part (4.4) 
and the rest 



1 



a = 



U{0)-^U{2kF) 



1 



U {2kF) 



= a. 



1 - 1 

AttVf Ns 



-91, 



where 



a. 



4nVF \^'~N^' 



Attvf 



(5.33) 



The condition for the tunneling exponent to be negative is more restrictive that it 
seemed to be: gi > Nsg2- It is not hard to see that the RG equation for now 
changes to 

? = -^ral-^^r,(0V (5.34) 



dl 

where Fi {I) is given by 



Integrating (5.34), we find 



Ti = 



'i-nvF Ns 
1 



NsQl W 



{\e\/Wr 



where 



N^ 



^as now electrons have spins, they are allowed to be at the same point in space and interact. 
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In particular, for Ng = 2, we get 



and conductance 

« = e"(' + 2£;'"F)'"'<W''»''"°--- 

where Go is the conductance for the free case. Thus the flow of the backscatter- 
ing amplitude results in a multipUcative log-renormaUzation of the transmission 
coefficient. One can check the first-order result is reproduced if expand the RG 
result to the first log. 

An interesting feature of this result is that it predicts three possible types of 
behavior of the conductance as function of energy. 

1. weak backscattering: 

a' > ^ 51 < g2/Ns. 

In this regime already the first-order correction corresponds to suppression of 
the conductance, which decreases monotonicaUy as the energy goes down. 

2. intermediate backscattering: 

a' >0 but a'^ < ^ < Qi < N^Qi- 

In this regime, the first-order correction enhances the transparency, but the RG 
result shows that the for e ^ 0, the transmission goes to zero. It means that 
at higher energies, when the RG has not set in yet, the conductance increases 
as the energy goes down, but at lower energies the conductance decreases. 
The dependence of G (e) on e is non-monotonic-there is a maximum at the 
intermediate energies. 

3. strong backscattering: 

In this regime, tunneling exponent a'^ is negative and the conductance in- 
creases as the energy goes down. 

5.4. Comparison of bulk and edge tunneling exponents 

Tunneling into the bulk of a ID system is described by the density of states ob- 
tained, e.g., in the DL solution of the Tomonaga-Luttinger model (no backscat- 
tering). The "bulk" turmeling exponent is equal to 

2„. 



4u 
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where 



u = a/1 + Igij-KVp. 

In this Section, we considered tunneling into the edge for weak interaction and 
found that the conductance scales with exponent 2a^ (5.33). To compare the two 
exponents, we need to expand the DL exponent for weak interaction 

2a = 



For gi = 0, the edge exponent is 

= 2^'^- 

We see that for weak coupling tunneling into the edge is stronger affected by the 
interaction than tunneling into the bulk: the former effect starts at the first order 
in the interaction whereas the latter starts at the second order. This difference has 
a simple physical reason which is general for all dimensions. In a translationally 
invariant system, the shape of the Green's function is modified in a non-trivial 
way only starting at the second order. For example, the imaginary part of the 
self-energy (decay of quasi-particles) occur only at the second order. The first- 
order corrections lead only to a shift in the chemical potential and, if the potential 
is of a finite-range, to a renormalization of the effective mass. If the translational 
invariance is broken, non-trivial changes in the Green's function occur already 
at the first order in the interaction. That tunneling into the bulk and edge are 
characterized by different exponents is also true in the strong-coupling case (cf . 
Sec.6). As the relation between the bulk and edge exponents is known for an 
arbitrary coupling, one can eliminate the unknown strength of interaction and 
express one exponent via the other. Knowing one exponent from the experiment, 
one can check if the observed value of the second exponent agrees with the data. 
This cross-check was cleverly used in the interpretation of the experiments on 
single- wall carbon nanotubes [71, 72]. 



6. Bosonization solution 

Bosonization procedure in described in a number of books and reviews [1]-[10]. 
Without repeating all standard manipulations, I will only emphasize the main 
steps in this Section, focusing on a couple of subtle points not usually discussed 
in the Uterature. A reader famiUar with bosonization may safely skip the first 
part of this Section, and go directly to Sees. 6.1.5 and 6.2.1, where tunneling 
exponents are calculated. Some technical details of the bosonization procedure 
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are presented in Appendix Appendix C. As in Sec. 3, = 1 in this Section, 
unless specified otherwise. 

6.1. Spinless fermions 

6.1.1. Bosonized Hamiltonian 

We start from a Hamiltonian of interacting fermions without spin 

p,k,q P^k,q 

The interacting part of the Hamiltonian can be re-written using chiral densities 
as 

Hint = 92^J2p+ P- + ;^ m ^+ P+ + P- P- ■ 

9 9 
The interacting part is in the already bosonized form. For a linearized dispersion 

= \k\ - kp, 

it is also possible to express the free part via densities. To check this, let's assume 
that Hq can indeed be written as 

Ho = jY^i ^+ (9) P+ (-9) + P- (9) P- (-9)] , 

9 

where Aq is some unknown function. Commuting with JJq, and making use 
of the anomalous commutator [p+ {q) , p+ {—q)] = qL/2-K, we obtain 

9' 

= \p+{q)-- 

TT 

On the other hand, the same commutator can be calculated directly in a model 
with the linearized spectrum, using only the fermionic anticommutation relations 
[69]. This gives 

[P+ (q) , Ho] = qp+ (g) . 
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Comparing the two results, we see that 



and thus 



-^0 = TT-^ ^ p+ {q) p+ {-q) + p- (q) p_ {-q) . 

Q 

Combining the free and interacting parts of the Hamiltonian, we obtain 

if = (l + |i) ^ {p+ (g) p+ {-q) + p_ (g) p_ (_g)}+7r^ff2 J2 P+ ('^) P- 

Q Q 

Notice that if only the 54— interaction is present, the system remains free but the 
Fermi velocity changes. 

It is convenient to expand the density operators over the normal modes 



p+{x) = ^^JiL(,,^e^«^+6te-^«x) 

q>0 



q<0 



One can readily make sure that density operators defined in this way reproduce 
the correct cormnutation relations, given that 
operators, the Hamiltonian reduces to 



5q^q' . In terms of these 



= ^ (1 + E ^{b% + blqb-q] + TTffa E « + ^h-q) . 



q>0 q 

Introducing new bosons via a BogoUubov transformation 

cl = cosh.6„h„ + smh.6„h^_, 



Jq"-q, 



cosh 6qb^_ + sinh 9qbq 



and choosing 9q so that the Hamiltonian becomes diagonal, i.e., 



tanh26i, 



52/27r 



we obtain 



" l + g4/2n' 
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where 



u\q\ , 



and 



31 



-,1/2 



(6.2) 



For the spinless case, backscattering can be absorbed into forward scattering. 
The resulting expression for the renormahzed velocity for the case 52 = 54 7^ 51 
is (cf. Appendix Appendix C.3 ) 



92-91 



2tt 



6.1.2. Bosonization offermionic operators 

The ^— operators of right/left movers can be represented as 



'21: a 



(6.3) 



where a is the ultraviolet cut-off in real space. Using the commutation relations 
for p±, one can show that the (anti) commutation relations 



{*± (a;),*^ (a;')} =5{x-x') 



are satisfied. 

The argument of the exponential can be re- written as two bosonic fields 



1 



V27ra 

f{x)T^ {x) ■ 



(6.4) 



Let's now introduce backscattering. Since for spinless fermions backscattering is just an ex- 
change process to forward scattering of fermions of opposite chirality (think of a diagram where 
you have right and left lines coming in and then interchange them at the exit), the only effect of 
backscattering is to replace 92 — > 92 — si . (cf. discussion in Sec. 4). Instead of (6.2), we then have 



/ 92 -91 
V 27r 



1/2 



(6.1) 



Now, consider a delta-function interaction, when gi =92 =94- Pauli principle says that we 
should get back to the Fermi gas in this case. However, u still differs from unity (Fermi velocity) 
and thus our result violates the Pauh principle. See Appendix Appendix C.3 for a resolution of the 
paradox. 
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Equating exponents in (6.3) and (6.4), we obtain 

Vtt [<p (x) - I? (x)] = 2tt ( dx'p+{x') 

J — OO 



(6.5) 



q>0 



2nL iq 

y/Tr[ip{x)+'d{x)] = 2n f dx'p_{x') (6.6) 

J — OO 

Solving for ip (x) and § {x) gives 

V{x) = -i Yl ^msgn9(6ge'«^-6te-"?-); (6.7) 

-00«J<CX) V^W^ 

-cx><g«x> V'^lll^ 

Using (6.7) and (6.8), one can prove that (p (x) and dx'& {x) satisfy canonical 
commutation relations between coordinate and momentum (cf. Appenidx Ap- 
pendix C.2.1) 

[ip{x),d^'^{x')]=i5{x-x'). (6.9) 

Using Eqs. (6.5) and (6.6), we obtain the density and current as the gradients of 
the bosonic fields 

/X pX 
dx' (p+ (a;') + p- {x')) = \ dx' p (x') — > 
-oo J —oo 

p {x) = -^dx(p; 

dx' {p+ {x') + p- {x')) = -V^ / dx'j {x') ^ 

-oo J — OO 

j {x) = ]=d^-& {x) . 

The continuity equation, 

dtp + dxj = 

relates the Heisenberg fields ip (x, t) and -d {x, t) 
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The current can be also found as 



j{x,t) = — -j=dtip{x,t) . 
We will use this relation later. Expressing H in the real space via p± 

H = TT j dx[pl + pV\+^ j dx [p\ + pV\ +92 j dxp+p. 
and using the relations 



P± = 



1^ 



we obtain a canonical form of H in terms of the bosonic fields 



H 



where 



dx 
dx 



92 +54 



dx {dx(p) 



92 - .94 

47r 



dx (dx^) 
(6.10) 



'1 + ^ 



1 + 



34+52 ■ 
2tt 



(6.11) 



(6.12) 



For 94 = 92 = g,we have 

u = y/TTgJ^; K = —2==. 

Notice that in this case uK = 1. This is important: in the next Section, we will 
see that this product renormalizes the Drude weight (and the persistent current). 
Neither of these quantities are supposed to be affected by the interactions, as the 
Galilean invariance remains intact. We see that it is indeed the case in our model. 

If backscattering is present (but 92 = 94), the parameters change to (cf. Ap- 
pendix Appendix C.3) 



u = x/l + 
K = 



92-91 



v/l + (.92-.gi)/7r 



(6.13) 
(6.14) 



As we have already seen in Sec. 3, a difference between 32 and 94 leads to the current-current 
interaction in the Hamiltonian. In the bosonized form, this interaction is the (92 — 94) {dx'&)^ term 
in the first line of Eq. (6.10). 
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Had we started with another microscopic model, e.g., with fermions on a lattice 
but away from half-filling, the effective low-energy theory would have also been 
described by Hamiltonian (6. 10) albeit with different-and, in general, unknown- 
parameters u and K. The term "Luttinger liquid" (LL) [70] refers to a universal 
Hamiltonian of type (6.10), which describes the low-energy properties of many 
seemingly different systems. In that sense, the LL is a ID analog of higher- 
dimensional Fermi Uquids, which describe the low-energy properties of a large 
class of fermionic systems, while encoding the quantitative differences in their 
high-energy properties by a relatively small set of parameters. 



6.1.3. Attractive interaction 

What happens for the case of an attractive interaction, 17 < 0? Formally, for 
g < — TT (or g2 — 91 < tt), in Eqs.(6.12,6.13) is negative, which seems to 
suggest some kind of an instability. Actually, this is not the case [59], as a ID 
system of spinless fermions does not have any phase transitions even at T — 0. 
All it means that the interacting system is a liquid rather than a gas, i.e., it does 
not require external pressure to mantain its volume. An equiUbrium value of the 
density is fixed by given ambient pressure. To see this, restore the Fermi velocity 
vp = TTn/m, where n is the density 

u^ = vl(l + ^]^(^)\^ (6.15) 

\ TTVF J \'m J TO 

and recall the thermodynamic relation 

w2 = m~^dP/dn, (6.16) 

where P is the pressure. Integrating (6.16) with the boundary condition P (n = 0) 
0, we obtain the constituency relation 

/7r\2n^ g 2 

P=\ — ) + • 

\mJ 3 2to 

For g < 0, there is a metastable region of negative pressure. This means that if 
the ambient pressure is equal to zero, the thermodynamically stable value of the 
density is given by the non-zero root of the equation P (n) = 0: 

n* = ^\g\m. 



The square of the sound velocity at this density is positive: 
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The Fermi velocity dXn = n* is 



v*p = -Kn* jm = \g\ 



73 



and parameter K 

U* TV 

is a universal number, independent of the interaction. 



6. 7.4. Lagrangian formulation 

In what follows, it will be more convenient to work in the Lagrangian rather the 
Hamiltonian formulation (and also in complex time). A switch from the Hamil- 
tonian to Lagrangian formulation is done via the usual canonical transformation 



(6.17) 



S = jdx jdt{qp-nip, q) ), 
where H is the Hamiltonian density defined such that 

H = j dxH 

and q and p are the canonical coordinate and momentum, correspondingly. Ac- 
cording to commutation relation (6.9), 



q = ipp = d^i^. 



(6.18) 



Performing a Wick rotation, t —ir, we reduce the quantum-mechanical prob- 
lem into a statistical-mechanics one with the partition function 



Z = j Dip j D-ffe-^", 



where the Euclidian action 

Se = I dr I dx 



In a Fourier-transformed form 

'1 u 



/I 1i 1 
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where k = {q,u)). If one needs only an average composed of fields of one type 
(ip or then the other field can be integrated out. This leads to two equivalent 
forms of the action 



-'- 2 I 2 

—uj + uq 
u 



/ dx 

2K 



dT 



u 



2 , 2 

—u> + uq 
u 



dT 



u 



(6.19a) 
(6.19b) 
(6.19c) 

(6.19d) 



In calculating certain correlation functions, e.g., the fermionic Green's func- 
tion, one also needs a cross-correlator . This one is computed by keeping 
both (f and d in the action. 

It is convenient to re-write the action in the matrix form 



where 



and the inverse matrix of propagators 



D- 



Inverting the matrix, we obtain 



q^uK 



D 



1 



iquj 
iqiv q^^ 



uK —iuj/q 

-iuj/q f 



The space-time propagators can be found by performing the Fourier transforms 
of D. For diagonal terms, one really does not need to do it, as it is obvious from 
(6. 19a) and (6. 19c) that these propagators just coincide with the Green's function 
of a 2D Laplace's equations. Recalling that the potential of a line charge is a 
log-function of the distance, we obtain 

K 



47r a;2 + [u |r| -|- a) 



(i?(^)i?(0)-i?2 (0)) 



■In- 
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where a is "lattice constant", z = {x,t) , and x'^ + v?t^ » c?. ^•^ These are 
the two correlation functions we will need the most. In addition, there is also an 
off-diagonal propagator 

S(^) = ((p (z) I? (0)-(^(0)^(0)) = (i?(z)^(0)- 19(0)^(0)) 



5 (2:) depends only on the ratio x/ut and thus does not change the power- 
counting. To see this, introduce polar coordinates q = k cos a/u, iv = sma,x = 
{z/u) cos (3, and t = zsm(3. Then 



1 j^J da (e"' - 1) tan a. 



(27r)' 7o 

The resulting integral is a function of only /? = tan~^ [x/ut). 
6.1.5. Correlation functions 

Now we can calculate various correlation functions, including the Green's func- 
tion. 

Non-time-ordered Green's function for right movers: 

G+{x,t) = -{T^^+{x,T)ijliO,0)) (6.20) 

= J_^jxBg«v^(¥>(l)-i?(l))g-iVi^(¥'(0)-i?(0))^^ 

where (1) = (x, r) and (0) = (.t = 0, t = 0) , and where is a bosonic time- 
ordering operator. I will use the well-known result, valid for an average of the 
product of the exponentials of gaussian fields (see the books by TsveUk [6] or 
Giamarchi [10] for a derivation) 

' " (6.21) 



^'A non-symmetric appearance of the cut-off with respect to time and space coordinates reflect an 
asymmetric way the sums over bosonic momenta and frequencies were cut. We adopted a standard 
procedure in which the sum of over q is regularized by cxp (— |fe| a) , whereas the frequency sum is 
unlimited. Other choices of regularization are possible. 

''^Surely, it is not a conventional definition of the Green's function, but it is easier to work with this 
one for now, and restore the fermionic Tr product at the end. 
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[Eq. (6.21) is essentially a field-theoretical analog of the probability theory result 
for the average of e''*''' , where 7 is a Gaussian random variable.] For example, in 
the average 

A-i = \pK, A2 = —\pK and 



xx^ + (u \t\ + of j 
Similarly, with the help of (6.21), Eq. (6.20) reduces to 

G+ {x,t) = J_e'^<'^(l)v'(0)-(p=(0)>^g7r{,?(l)^(0)-,?^(0)),g-2{¥>(l)!?(0)-v(0)tf(0)>, 
_ ^7i-$(a;,r)^7re(z,T)^-2H(a;,r) (6.22) 

27ra 




2 \ 



+ (u |t| + a) 



e ■ 



where {■ ■ ■)t stands for a time-ordered product and where it was used that in a 

translationally invariant and equilibrium system (95^ (0)) = {ip^ (1)) (same for 
§) . Function / (x/mt) is a phase factor which does not effect the power-counting. 

Bulk tunneling DoS For a; = 0, 



By power-counting. 



G(0,t) a T-^^ 



u (e) a |£| 2 i = |£| . (6.23) 



This is an analog of the DL result for the spinless case. 

Edge tunneling DoS In a tunneling experiment, one effectively measures the 
local DoS at the sample's surface. In a correlated electron system, the boundary 
condition affects the wavefunction over a long (exceeding the electron wave- 
length) distance from the surface. Therefore, the surface DoS differs signifi- 
cantly from the "bulk" one. If a tunneling barrier is high, then-to leading order 
in transmission- the DoS can be found via imposing a hard-wall boundary con- 
dition. The presence of the surface (boundary) can be taken into account by 
imposing the boundary conditions on the number current 

j (a; = 0, r) = - ^drf = 0. (6.24) 
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at a; = 0. This means that is pinned at the boundary, i.e., it takes some time- 
independent value. In the gradient-invariant theory, we can always choose this 
constant to be zero. Thus, 

^(0,r) = 0. 

This suggests that the local correlator $(0, t) = 0, and the long-time behavior 
of the Green's function in Eq. (6.22) is determined only by the correlator of the 
1? fields. Had the boundary not affected this correlator, we would have arrived at 

G (a; = a;' = 0, r) (X exp (ttO {x = 0, r)) ex -^-^J2K' (wrong) 

But then we have a problem, as Eq. (wrong) does not reproduce the free-fermion 
behavior for iiT = 1. Consequently, the DoS at the edge (e) oc |e| would 
have not reproduced the free behavior either. What went wrong is that we pinned 
one field but forgot the other one is canonical conjugate to the first one. By 
the uncertainty principle, fixing the "coordinate" ((^) increases the uncertainty 
in the "momentum" (t?)-and vice versa. Thus, fluctuations of 'd fields should 
increase. A rigorous solution to this problem is to change the fermionic basis 
from the plane waves to the solutions of the Schrodinger equation with the hard- 
wall boundary condition and to bosonize in this basis. This was done by Eggert 
and Affleck [75] and Fabrizio and Gogolin [76] , [9]. Here I will give an heuristic 
argument based on a simple image construction, which leads to the same result. 
Eq. (6.24) translates into the boundary conditions for the bosonic propagators: 

^e{x,x',T)=Q; a^,^/ee(a;,a;',a;) = 0, (6.25) 

for X, x' = 0, where subindex e denotes the correlators in a semi-infinite system. 
Since $e and 9e satisfy the Laplace's equation, we can view these propagators 
as potentials produced by some fictitious charges. Then, <i?e and 6e can be con- 
structed from the propagators of an infinite sample by the method of images: 

^e{x,x',T) = ^{x - x' ,t) - ^{x + x',t); 

Qe{x,x',T) = Q{x-x',t) + Q{x + x',t). 

For X = x', 

$e(0,0,r) = 0; ee(0,0,T) = 2e(0,T). (6.26) 

Hence, pinning the (p field enhances the rms fluctuations of the i? field by a factor 
of two. This leads us to 



G+(0,0,r) cx exp(7r$e(0,0,r))exp(7ree(0,0,T)) 

' 2tt a 
2^ ""k 



exp (27ree (0, r)) oc exp ( In ) oc |r| . 
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Consequently, the DOS becomes 

Ve (e) oc |£|^"'"^ . (6.27) 

This result by Kane and Fisher [68] initiated the new (and still continuing) surge 
of interest to ID systems (in terms of the impurity scattering time, this result was 
obtained earher in Refs. [73, 74]). For tunneling from a contact with energy- 
independent DoS ("Fermi Uquid") into a ID system, the tunneling conductance 
scales as (e) 

Q{£) oc fe(e) oc \£\^ ~^ . 

Now we see that the free-fermion behavior is correctly reproduced for = 1. 

Expanding the tunneling exponent — 1 with parameter K from Eq. (6.14) 
for the weak-coupUng case gives 

j^-i _ ^ _ 92 -gi 
2'kvf 

This is the same result as the weak-coupling tunneling exponent (5.30) obtained 
in Sec. 5 via the scattering theory for interacting fermions. 

Where do the "bulk" and "edge" forms of DoS match? Consider an object 
G{x = x',£). At the boundary, the DoS is of the "edge" form (6.27). Far 
away from the boundary, the Green's function does not depend on x and v{e) 
acquires a "bulk" form (6.23). As a function of x, G{x = x' ,e) varies on the 
scale ~ u/\£\ and the crossover between two Umiting forms of v occurs on this 
scale. Choosing the energy in a tunneling experiment, i.e., temperature or bias- 
whichever is larger, determines how far from the boundary one should go in order 
to see a change in the scaUng behavior. 

6.2. Fermions with spin 

For fermions with spin, each component of the fermionic operator is bosonized 
separately 

V'±,<T = — exp [±iy/T:{ipa- T i?<t)] • 
V27ra 

Index (T of the bosonic field does not mean that bosons acquired spin. We simply 
have more bosonic fields. Charge and spin densities and currents are related to 
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the derivatives of the bosonic fields 

P±,<T = TT-l^iv'a^K) P<y^ P+,<T+ P-,<T = A=v'a'^ 

3a = P+,a - P-,a= -^K'^ 
1 

Pc 



ps = PT - Pi = ^ K - <^1) = <P'. 



where ' denotes dx and where the charge and spin bosons are defined as 

— V5 — — (o-26) 

I assume that the interaction is spin-invariant, i.e., couplings of tt and ti fermions 
are the same. Substituting the relations between charge- and spin-densities into 
the Hamiltonian, one arrives at the familiar bosonized Hamiltonian which con- 
sists of totally independent charge and spin parts 

H = Hr + Hg, : 



+ -^^ [ dxcos(V8Tr(t)a) . (6.29) 
(27ra) J ^ 

Parameters of the Gaussian parts are related to the microscopic parameters of the 
original Hamiltonian 

1/2 / 1 , /O^ \ 1/2 



I giy/Vi I 4ff2-gi Y^\ ^ _ 1^ i + gi/27r y 

V+2n) 27r J ' ' [l + {Ag2 - gi) /2n J 

m 



277 y y ' vi- 51/277 
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Notice that <l for g\ < 2g2 ("repulsion") and Kc> 1 for gi > 2g2 ("attrac- 
tion"). The boundaries for "repulsive" and "attractive" behaviors coincide with 
those obtained when studying tunneling of interacting electrons. The velocity of 
the charge part for gi =0 coincides with that found in the DL solution (Sec. 3) 



Scaling dimension of the backscattering term in the spin part can be read off 
from the correlation function 



If we allowed for different couphng constants between electrons of different spin 
orientations, then the coefficient in front of the cos term would have been gi±. 
For Ks > 1, the operator scales down to zero as a ^ 0, whereas for < 1, it 
blows up signaling an instabihty: a spin-gap phase. 

The RG-flow of the spin-partis described by the Berezinskii-KosterUtz-Thouless 
phase diagram. The fixed-point value of = for K* > 1. In the weak cou- 
pling hmit, the RG reduces to a single equation for gi , which we have derived in 
the fermionic language in Sec .4 

6.2.1. Tunneling density of states 

The procedure of finding the scahng behavior for the DoS reduces to a simple 
recipe. 

• Take the free Green's function and split it formally into the spin and charge 
parts 

1 _ 1 1 

• In an interacting system, the exponent of 1/2 in the charge part is replaced by 

(Kc + K-'^) /4 and in the spin-part by (Ks + K''^) /4. If the spin-rotational 
invariance is preserved, then the spin exponent remains equal to 1/2. 

• For a; = 0, 

G (t) oc 

and the DoS behaves as 
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Comparing this result for gi = Q with that by DL (Sec. 3), we see that the 
bosonization solution gives the same result as the fermionic one. 
• For tunneling into the edge, remove K^, which comes from the correlator 
{ipip) pinned by the boundary, and multiply K~^, which comes from {'&'d), by a 
factor of 2. This gives 

Ge it) OC 

and 

g OC I/e(£) OC |£|'(^="'"^) . 
Expanding Kc back in the interaction 

^ ^ / l + ffi/27r y^\^ g2-(l/2)ffi 
\l + {4g2-gi)/2nj ~ tt 

we obtain the weak-coupUng limit for the tunneling exponent 

(l/2)(l/i..-l«^^^). 

This coincides with the result obtained in the fermionic language (Sec.5). What 
was missed in a bosonization solution is a multiplicative log-renormalization, 
present in Eq. (5.35). This is because we evaluated G at the fixed point, where 
gl = 0, rather then derived an independent RG equation for the flow of the 
conductance. This procedure should bring in the log-factors (cf. Ref.[74] where 
these factors were obtained for the impurity scattering time). 



7. Transport in quantum wires 

7. 1 . Conductivity and conductance 
7.1.1. Galilean invariance 

Interactions between electrons cannot change the response to an electric field in 
a Gahlean-invariant system-the electric field couples only to the center-of-mass 
whose motion is not affected by the inter-electron interaction. This property is 
reproduced by the bosonized theory provided that the product uK = 1 (= f^? in 
dimensional form.) To see this, combine the Heisenberg equation of motion for 
density p (spinless fermions) with the continuity equation: 

dtp = i[H,p] = -d,j. (7.1) 
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Calculating the commutator in Eq. (7.1) with the help of Eq. (6.9), we identify 
the current operator as 



The current is not affected by the interaction as long as uK = 1. 



7.1.2. Kubo formula for conductivity 

The Kubo formula relates the conductivity, a response function to an electric field 
at finite w and q, to the current-current correlation function 



o- (w, q) = 

ILO 



— + {JJ) 
TT 



(7.2) 



where it was used that n = kp/n and kp/m = vp = lin our units. 
Electric current for electrons (e > 0) 



In complex time. 



e 

I 



i-dl) {m),,r 



e2 e2 



TT + U q TT 



(7.3) 



The first term in (7.3) cancels the diamagnetic response in (7.2). Continuing 
analytically to real frequencies, we find 



1 



n iu}-(uj + i5) + u^q^ 



(7.4) 
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Consequently, the dissipative conductivity is equal to 



1 

Recr {w,q) = wlm- 



TT — u^g^ + isgnuS 

— [5{io- uq) + (5 (w + uq)] . (7.5) 



7.1.3. Drude conductivity 

In a macroscopic system, one is accustomed to take the limit g — > first: this cor- 
responds to applying a spatially uniform but time-dependent electric field [61]. 
(For the lack of a better name, I will refer to the conductivity obtained in this way 
as to the Drude conductivity). The Drude conductivity in our case is the same as 
for the Fermi gas as the charge velocity drops out from the result 

Recr (w, 0) = 6^5 (w) 

or, restoring the units, 

2 

Re(7(a;,0) = ^^5{u). 

All it means that when a static electric field is appUed to a continuous system of 
either free or interacting electrons, the center-of-mass moves with an acceleration 
and there is no linear response, as there is no "friction" that can balance the 
electric force. 

For electrons with spins, the electrical current is related only to the charge 
component of the {}— field: 

Jc = -ejc = C'^dx'dc, 

where again UcKc = 1. Because of the a/2 factor in the current, the conductivity 
is by a factor of two different from that in the spinless case 

Rea(w,0) = ^^-^S{u;). 



(Notice, however, that at fixed density v i? is by a factor of 2 smaller for electrons 
with spin, so that the conductivity is the same.) 
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7.1.4. Landauer conductivity 

Let's consider now the opposite order of limits, corresponding to a situation when 
a static electric field is applied over a part of the infinite wire. (Again, for the lack 
of a better name, I will refer to this conductivity as to the Landauer conductivity.) 
The electric field might as well be non-uniform; the only constraint we are going 
to use is that the integral 



/ 



dxE [x) , 



equal to the applied voltage, is finite. The induced current (which in ID coincides 
with the current density) is given by 



J{t,x) = j dx' j dt'a{t-t';x,x')E{t',x') 



dx' I '■^e-'^'a{uj;x,x')E{oj,x'). 

J 277 

In linear response, the conductivity is defined in the absence of the field. As such, 
it is still a property of a translationally invariant system and depends thus only on 
X — x'. This allows one to switch to Fourier transforms 

J{t,x) = Jdx'J^J ^e<^--')e-'^*aico,q)E{w,x'). (7.6) 

Now use the fact that the appUed field is static: E {x, w) = 2^5 (w) Eq (x) (upon 
which the f-dependence of the current disappears, as it should be in the steady 
state) 

J{x) = jdx' j ^e'<--^')a (0, q) Eo (x') . (7.7) 

From (7.5), 

a{0,q) = -e^S{q) = Ke^S{q), (7.8) 
u 

where uK = 1 was used again. Substituting (7.8) into (7.7), we see that the x— 
dependence of the current also disappears 

Conductance Q = J/V is given by 



2n ' 
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or, restoring the units, 




(7.9) 



For electrons with spin, a similar consideration gives 




(7.10) 



We see that the conductance is renormalized by the interactions from it universal 
value given by the Landauer formula for an ideal wire [68]. 

7.2. Dissipation in a contactless measurement 

What kind of an experiment Eqs. (7.9) and (7.10) correspond to? 

Suppose that we connect a wire of length L to an external resistor and place 
the whole circuit into a resonator [79]. Now, we apply an ac electric field E (x,t) 
of frequency and parallel to a segment of the wire of length Le <^ L, and 
measure the losses in the resonator. The external resistor takes care of energy 
dissipation: as the wire is ballistic (also in a sense that electrons travel through 
the wire without emitting phonons), the Joule heat can be generated only outside 
the wire. Dissipated energy, averaged over many periods of the field, is given by 



For a monochromatic field, E (x, t) = Eq (x) coswoi and after averaging over 
many periods of oscillations, we obtain 



where u is the velocity of the charge mode in the wire. Because the wavelength 
of the charge excitations at frequency cjo-^coustic plasmons- is much shorter 
than the distance to contacts (L), the conductivity is essentially the same as for 
an infinite wire and depends only on a; — a;'. Performing partial Fourier transform 
in Eq. (7.5), we find 





Now, choose the frequency in such a way that 



Le <^ — <^ L 



(7.11) 



Recr (w, x) = - — cos (ujx/u) = —K cos {uix/u) , 



a. 12) 
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so that 

(5 = — j '^^ j '^2;'cos[a;o {x — x') /u]Eq{x) Eq{x') . 

On the other hand, because \x,x'\ < Le <^ u/uq, the cosine can be replaced by 
unity, and 

Q = ~KV^ = -GV\ 

or 




Therefore, dissipation in a contactless measurements under the conditions spec- 
ified by Eq. (7.11) corresponds to a renormalized conductance. To the best of 
my knowledge, this experiment has not been performed. A typical (two-probe) 
transport measurement is done by applying the current and measuring the voltage 
drop between the reservoirs. In this case, the measured conductance does not cor- 
respond to Eqs.(7.9,7.10) but is rather given simply by e^//i per spin projection- 
regardless of the interaction in the wire [80], [81], [82]. This result is discussed in 
the next Section. 

7.3. Conductance of a wire attached to reservoirs 



The reason why the two-terminal conductance is not renormalized by the inter- 
actions within the wire is very simple. For the Fermi-gas case, the conductance 
of e^//i per channel is actually not the conductance of wire itself-a disorder-free 
wire by itself does not provide any resistance to the current. In a four-probe mea- 
surement, when the voltage and current are apphed to and measured in different 
contacts, the conductance of a disorder-free wire is, in fact, infinite. However, in 
a two-probe measurement, the voltage and current contacts are the same. Finite 
resistance comes only from scattering of electrons from the boundary regions, 
connecting wide reservoirs to the narrow wire [83, 84], as shown in Fig. 15a). 
The universal value of e^/h is approached in the limit of an adiabatic (smooth on 
the scale of the electron wavelength) connection between the reservoirs and the 
wire [85] As the resistance comes from the regions exterior to the wire, the 
interaction within the wire is not going to modify the e'^/h— result. Another way 
to think about it is to notice that the renormalized conductance (7.9,7.10) can be 

'^Accidentally, the actual constraint on the adiabaticity of the connection is rather soft-it is enough 
to require the radius of curvature of the transition region be just comparable to, rather than much larger 
than, the electron wavelength [85]. 
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Fig. 15. a) A Luttinger-liquid (LL) wire attached to Fermi-liquid (FL) reservoirs, b) same for a single 
impurity within the wire. 

interpreted as a manifestation of the fractional charge e* = \fK{^JKc), associ- 
ated with the excitations in a ID system. However, the current coming from, e.g., 
the left reservoir is carried by integer charges, and as all these charges get even- 
tually transmitted through the wire, the current collected in the right reservoir 
is carried again by integer charges. Fractional charges is a transient phenomena 
which, in principle, can be observed in an ac conductance or noise measurements 
but not in a dc experiment. In the rest of this Section, these arguments will be 
substantiated with some simple calculations. 

7.3.1. Inhomogeneous Luttinger-liquid model 

An actual system consists of two Fermi-liquid reservoirs connected via a Luttinger- 
liquid (LL) wire and , due to the presence of the reservoirs, is not one-dimensional. 
In the inhomogeneous Luttinger-liquid model, the actual system is replaced by an 
effective ID system, which is an infinite LL with a position-dependent interac- 
tion parameter K{x) (cf. Fig. 16). The actual reservoirs are higher (£> = 2 or 
3) systems, where the effect of the interaction can be disregarded. Consequently, 
the reservoirs are modelled by one-dimensional free conductors with A'l = 1 . 
In between, K (x) goes through some variation. Similarly, the charge velocity 
is equal to the Fermi one in the reservoirs and varies with x in the middle of the 
system. The potential difference applied to the system produces some distribu- 
tion of the electric field along the wire. The shape of this distribution is irrelevant 
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K(x) 



K=1 

L 



V(X) 



E(x) 




Fig. 16. Inhomogeneous Luttinger-liquid model. 



in the dc linear response. 

7.3.2. Elastic-string analogy 

The (real-time) bosonic action for a spinless LL is 



1 1 



1 



2 J K(x) { u[x) 



(7.13) 



The density of the electrons (minus the background density) and the (number) 
current are given by 

p = da^^pl^fii, j = ~dt(p/Vn. (7.14) 
The interaction with an external electromagnetic field is described by 



so that the equation of motion for is 



(7.15) 



where E = —OxAq + dtAi is the electric field. We assume that the electric field 
is switched on at f = 0, so that E{x,t) = for t < and E{x,t) ~ E{x) 
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Fig. 17. a) Solution of the wave equation in the homogeneous case for t = bL/u. b) Schematic 
solution in the inhomogeneous case for t^L/u. 



for t > 0. The problem reduces now to determining the profile of an infinite 
elastic string under the external force. In this language, tf{x, t) is the transverse 
displacement of the string at point x and at time t, while the number current 
j — —dt^p/ \pn is proportional to the transverse velocity of the string. 

To develop some intuition into the solution of Eq. (7.16), we first solve it in the 
homogeneous case, when K =const, u =const, and E(x) =const for |a;| < i/2, 
and is equal to zero otherwise. In this case, the solution of Eq. (7.16)fort > Lju, 
is 

f t-^;^for|a;|<L/2; 

, KeV t-\x\luioTLl2<\x\<ut~Ll2 

w[x, t) — ■= X < I I ' 'L, 9 /-II- / 

' 2V^ I ^(t- N^^2 ^orMt-L/2< < ui + L/2 

I " 0, for |x| > ut + L/2, 

where V ~ EL is the total voltage drop. This solution is depicted in Fig. 17a. 

The profile of the string consists of two segments (I and II in Fig. 17a) whose 
widths, equal to [ut — L), grow with time, and of three segments (III, IV and V in 
Fig. 17a) whose widths are constant in time and equal to i. In segments I and II, 
the profile of the string (p{x, t) is linear in x, and therefore, being the solution of 
the wave equation, also in t; in segments III-V, the profile is parabolic. Outside 
segments IV and V, the string is not perturbed yet, and (p{x,t) = 0. As time 
goes on, the larger and larger part of the profile becomes linear. For late times, 
the pulse produced by the force spreads outwards with velocity u, involving the 
yet unperturbed parts of the string in motion; simultaneously, in all but narrow 
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segments in the middle and at the leading edges of the pulse, the string moves 
upwards with the t- and a;-independent "velocity" dtf = KeV /2^Jtt. In terms 
of the original transport problem, it means that the charge current J = — ej 
is constant outside the wire (but not too close to the edges of the regions of 
where the electron density is not yet perturbed by the electric field) and given by 
J = Ke^V/h. Therefore, the conductance (per spin orientation) isQ = Ke^ /h. 

We now turn to the inhomogeneous case. As in the previous case, the profile 
consists of several characteristic segments (cf. Fig. 17b). In segments III-V, the 
profile is affected by the inhomogeneities in K{x), u{x), and E(x) and depends 
on the particular choice of the .x-dependences in all these quantities. In segments 
I and II however, the profile, being the solution of the free wave equation, is 
again linear in x (and in t). Requiring the slopes of the string be equal and 
opposite in segments 1 and II (which is consistent with the condition of the current 
conservation), the solution in these regions can be written as <^(a;, t) = A{t — 
\x\/ui). The constant A can be found by integrating Eq. (7.16) between two 
symmetric points ±a , chosen outside the wire 

-j_^dxd.{-d.^)=-^j_JxEix) = -j=. a. 11) 

Outside the wire, K{x) — Ki^ and u{x) = ul, thus A = Ki^eV/2^^h. Calcu- 
lating the current, we get Q = Ki^e^/h and, recalling that i^L = 1, we finally 
arrive at C/ = e'^/h. Thus, the conductance is not renormalized by the interactions 
in the wire. 

7.3.3. Kubo formula for a wire attached to reservoirs 

The Kubo formula for a translationally non-invariant system can be written as 

u{u);x,x') = 5 [x — x') 

inui 

+ [ diT-T')e'^--{TrdrV{x,T)dr'^{x',T'))\u 

The diamagnetic contribution is cancelled by a delta-function term, which is ob- 
tained when integrating by parts in the time-ordered product [86, 10]. Having 
this in mind, I will re-write the conductivity via the Fourier transform of the tpf- 
correlator without the T— product 

e2 

a (w; X, x') = i — {x, x') \iu;^-,oj+i6- 

TTU) 

For a translationally invariant case, this reduces back to Eq. (7.4). Now, K (x) 
and u {x) depend on position. The propagator of the ip fields satisfy the wave 
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equation (or a Laplace's equation as we are dealing with the imaginary time) 

^^^{x,x')=5{x-x'). (7.18) 



a. ( 4Sa 



_u{x)K{x) \K{x) 

In a model of step-like variation of K (x) and u {x) (K = Kw and u = uw 
within the wire and K = Kl = 1 and u = ul = I outside the wire), Eq. 
(7.18) is complemented by the following boundary conditions: 1) {x, x') 
is continuous at a; = ±L/2, 2) u {x) K {x) d^^uj^ {^i x') is continuous at a; = 
±L/ 2; but 3) undergoes a jump of unit height at x = x' . Solution of this problem 
is totally equivalent to finding a potential of a point charge located somewhere 
in a sandwich-like system, consisting of three insulators with different dielectric 
constants. Two of these layers are semi-infinite and the third one (in the middle) 
is of finite thickness. "Potential" can be found in a general form for 

arbitrary a;, a;'. In the expression for the current 

J{t,x) = J J '^'^ {'^'^x^x') E {uj,x') , 

x' is within the wire; hence we need to know (x, a;') only for —L/2 <x'< 
L/2. In a steady-state regime, one is free to measure the current through any 
cross-section; let's choose x also within the wire. As we are interested in the 
limit a; — > 0, when the plasmon wavelength is larger than the wire length, we can 
put X = x'. In the interval —L/2 < x = x' < L/2 the solution of the Laplace's 
equation is 

$ (xx)- + ^^-e-^/^" + ^+^- cosh(2.x/LJ 

where La; = uw / l^m | , uw is the charge velocity within the wire, and 

Letting uj^ in (7.19) to zero (and thus L^ to oo), we find that 

2iuj 2iuj 

as by our assumption Kl = 1. This result is true for any x, x' within the wire for 
w ^ 

(a;, a;') = — , for - L/2 < x,x' < L/2. 
2iuj 

The Luttinger-liquid parameters of the wire drop out from the answer. The dc 
conductivity reduces to its free value 
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and, consequently, the conductance 



is not renormalized by the interaction. The same consideration for electrons with 
spins gives 

Q = ^ (7.20) 

h 



7.3.4. Experiment 

Most of the experiments on quantum wires indeed show that the conductance is 
quantizated in units of 2e^//i at relatively high temperatures. At lower 

temperatures, the conductance decreases beyond the universal value and also the 
quantization plateaux exhibit some structure as a function of the gate voltage 
[90]. This can be interpreted as the effect of residual disorder: as was discussed 
in Sees. 5,6 transmission decreases at lower energy scales. An effect of single 
impurity in a quantum wire will be largely insensitive to the presence of reser- 
voirs: as long as the transmission coefficient for an impurity is much smaller than 
one, the largest voltage drops occurs near the impurity rather than at the contacts 
to the wire. One can show that the scaling of the conductance with energy is de- 
termined by the interaction parameter K inside the wire [91], in a contrast to the 
disorder-free case when only K outside the wire matters. Also, the mesoscopic 
conductance fluctuations increase as the temperature goes down (the theory pre- 
dicts that this effect is enhanced by the interaction [92]). As one is dealing here 
with a crossover regime from scattering at a single impurity to that at many impu- 
rities, a quantitative analysis of the temperature dependences is difficult; another 
compUcation arises from the finite length of the wire which cuts off scalings with 
temperature and voltage. In addition, at higher temperatures the first quantization 
plateau exhibits a well-defined step at about 0.7 x 2e'^/h [93, 94, 95, 96, 97, 98]. 
This "0.7" feature is not likely to result from spurious impurity scattering but 
rather reveals some interesting physics beyond what has been discussed so far in 
this review. Although the "0.7" feature deserves a review on its own, I will come 
back to this subject briefly in the next Section. 



However, it has been observed recently that the conductance of carbon nanotubes is quantized in 
units of /h -as opposed to ie^ /h, predicted by the non-interacting theory for this case [87]. 

'^A special case of a non-universal conductance quantization is very long wires grown by cleaved- 
edge overgrowth technique [88] can be attributed to a non-trivial coupling between the wire and 2D 
reservoirs [89], characteristic for these systems. 
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7.4. Spin component of the conductance 

As we have shown in the previous sections, the Luttinger-liquid models predicts 
that conductance of a disorder- free wire is given by /h per channel at any tem- 
perature. Also, thanks to spin-charge separation, spin degrees of freedom do not 
play an essential role in charge transport except for giving an overall factor of 
two to the conductance. These two results hold as long as the Luttinger-liquid 
model is a good description for interacting electrons in the wire. When does this 
model break down? If the interaction is strong, electrons form almost a periodic 
ID structure: quasi- Wigner crystal. The exchange energy of almost localized 
electrons is exponentially small and, correspondingly, the spin velocity is small 
too: Us <C Uc- The Luttinger-liquid model should work for energies (temper- 
atures) much smaller than the smallest of the two (spin and charge) bandwidth 
T <C Uskp ^ Uckp^ when both spin- and charge degrees of freedom are coher- 
ent. The spin- incoherent regime, i.e., Uskp <IC T <IC Uckp, has attracted con- 
siderable interest recently [19, 20, 21], and was shown to spoil the conductance 
quantization in integer multiples of 2e^/h [19] at temperatures larger than the 
spin bandwidth (Uskp). In what follows, I present a short summary of Ref.[19]. 

In a quasi- Wigner-crystal regime, a reasonable starting point for describing 
the spin sector is the Heisenberg model 



where spins are localized at "lattice sites" corresponding to positions of electrons. 
Because the Lieb-Mattis theorem [46] forbids ferromagnetic ordering in ID, the 
sign of the exchange interaction must be antiferromagnetic: Jex > 0. Assuming 
that electrons are well localized at distances a = l/n from each other, J^x can be 
estimated in the WKB approximation: J^^'^ Ep exp (— c/ y/aBfi) , where c ~ 1 
and ttB is the Bohr radius. A spin- 1/2 chain is then mapped onto a Hubbard 1/2- 
fiUed model of spinless fermions via the Jordan- Wigner transformation 




a]ai - 1/2; 



Sx (l) + iSy (/) = a] exp ( in ^ 




with the result 




The spinless Hubbard model can be bosonized 




(7.21) 
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A comparison to the Bethe Ansatz solution of the Hubbard model at half-filling 
enables one to identify the parameters of the spinless LL with the microscopic 
parameters of the spin- 1/2 chain. In particular, for an isotropic spin chain {Jx = 

Jy = Jz) 

u = I Jexa; K = 1/2. (7.22) 

Comparing the spin-part of the Hamiltonian of the original LL model, Eq. (6.29) 
with that of the spinless LL, Eq. (7.21), one notices that they the same upon the 
following mapping 

</) = -^<^«; e = V20s; ^ = UsKs = 2uK, (7.23) 

or 

Us = u^ I Jexa; Ks = 2K=l. (7.24) 

As Jex is exponentially small, so is the spin bandwidth. Therefore, the Luttinger- 
liquid description is valid only at very low temperatures. 

A translationally invariant LL still possesses spin-charge separation. However, 
this is no longer true for a wire connected to non-interacting leads. To understand 
this point, let's come back to the inhomogeneous LL model (cf. Sec. 7.3.1), 
where the electron density changes from a higher value in the leads to a lower 
value within the wire. Because Jex depends on the local density, it is modulated 
along the wire, and its minimum value is at the middle of the wire. In the leads, 
we have a non-interacting system, where ,/ex ^ Ep ^ T. However, in the middle 
of the wire spins are incoherent, if J™'" <C T. Thus a spin part of the electron 
incoming from the lead at energy T above the Fermi level cannot propagate freely 
through the wire because the spin band narrows down: it works as if there is a 
barrier for spin excitations in the wire. Although charge plasmons propagate 
freely, backscattering of spin plasmons leads to additional dissipation, and thus 
to additional resistance. The total resistance of the wire consists now of two parts 

R = i?c + Rs ■ 

The charge part, i?c, is due to propagation of charge plasmons. Since the charge 
part is still described by the LL model, our previous result for universal con- 
ductance, Eq. (7.20), still holds and Re = G''^ = h/2e^. For T < J™'", only 
athermal spin plasmons, with energies exceeding the width of the spin band, con- 
tribute to Rs - The number of such plasmons is exponentially small, hence 

Rs a exp {-J^'^/T) , 

and total conductance Q = (i?c + -Rs)~^ is exponentially close to 2e^/ h. At high 
temperatures (T J^'"), almost all spin plasmons are reflected by the wire. 
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Then Rg ^ Rc h/e^, and the conductance differ substantially from its univer- 
sal value. This qualitative picture can be confirmed in a particular simple case of 
the XY— model for the spin-chain. In this case, Rs can be calculated explicitly 
with the result [19] 

R ^ ^ 

' 2e2exp(J™/r) + l 

and, consequently, the conductance is equal to 

l + [exp(J™m + !]"'■ 

For T ^ 0,g approaches the universal value of 2e'^/h. For T > Je™'°,C? ap- 
proaches another T -independent Umit, equal to {2/'i)2e^/h. The actual number 
in the high-temperature limit of the conductance is model-dependent (it is dif- 
ferent, for example, for an isotropic spin-chain), but the main result, i.e., the 
non-universality of conductance quantization at higher temperatures, survives. 

As it was mentioned in Sec. 7.3.4, the experiment shows that there is a shoul- 
der in the conductance preceding the first quantization plateau at a fractional 
value of about 0.7x2e'^/h. Surprisingly, this "0.7 feature" is more pronounced 
at higher temperatures, and the T- dependence of this feature was reported to be 
of an activated type [95]. The magnetic field transforms the "0.7 feature" into a 
fully developed quantization plateau at e^/h, which is to be expected in a fully 
polarized, and thus spinless, regime. The sensitivity to the magnetic field hints 
at the spin origin of the effect, and a significant theoretical effort was invested 
in understanding how spins can explain the observed phenomena. Although the 
effect, described in this Section, does have all qualitative characteristics of the 
observed "0.7 feature", it is not clear at the moment whether this feature indeed 
corresponds to the spin-incoherent regime. Other explanations of the "0.7 fea- 
ture" have been suggested (most prominently, the Kondo physics is believed to 
be involved [99, 98]), but a further discussion of this point goes beyond the scope 
of these notes. 

7.5. Thermal conductance: Fabry-Perrot resonances of plasmons 

There is an important difference between charge and thermal (electronic) con- 
ductances [22]. As we have just shown, the charge conductance is equal to e^/Zi 
regardless of interaction in the wire. This means that the transmission coefficient 
of electrons is equal to unity. The effect of the temperature on the charge con- 
ductance is the same as for a non-interacting, perfectly transmitting wire: at finite 
temperature, not only the lowest but also higher subbands of transverse quanti- 
zation are populated, and the quantization plateaux are smeared. However, this 
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effect is exponentially small for temperatures smaller than either the Fermi en- 
ergy, Ep , or the difference between the Fermi energy and the threshold of the 
next subband of transverse quantization, A; whichever is smaller. 

Thermal current is carried not by electrons but bosonic excitations: acous- 
tic plasmons. In contrast to electrons, plasmons get reflected at the boundary 
between the reservoirs and the wire due to the mismatch of charge velocities 
(this reflection happens even for an adiabatically smooth transition). From the 
plasmon's point-of-view, a wire coupled to reservoirs represents a Fabry-Perrot 
interferometer Interference of plasmon waves scattered off the opposite ends of 
the wire results in an oscillatory dependence of the transmission coefficient on 
the frequency with a period given by the travel time of plasmons through the wire 

2'KU)L = L/uw- 

As long as Af <C L, this period is long: wl <C Ep ■ The difference between 
charge and heat transport is that the chemical potential of plasmons is equal to 
zero, and thus the characteristic scale for frequency is set by T. Therefore, the 
thermal conductance varies with the temperature on a scale T ~ ujl- 

Suppose that a small temperature difference 5T is maintained between the 
reservoirs, connected by a quantum wire. As the Hamiltonian of an interacting 
system is diagonalized of terms of plasmons, plasmons modes are decoupled and 
contribute to the energy flux independently. Then the thermally averaged energy 
current, i.e., the thermal current can be found via a Landauer-like argument 

where nL,R i^^/T) are the Bose distribution functions in the reservoirs. Expand- 
ing in small 6T, we obtain the thermal conductance 

For a free system, \t (w)]^ = 1 and Gt = nT/6. The charge and thermal con- 
ductances of a free system obey the Wiedemann-Franz law 



Qt 



Lo = —, (7.26) 



where Lq is the Lorentz number. This means that charge and energy are carried 
by the same excitations. This is not so in a Luttinger liquid. 

For an interacting system, relation (7.26) holds in the limit of T ^ 0. The 
characteristic scale for frequencies in integral (7.25) is determined by T. For 
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T <^ u)L, one can substitute a; = into \t ■ Regardless of the interaction 
strength, |t(0)|^ ~ 1 : a Fabry-Perrot interferometer becomes transparent in the 
long wavelength limit. For T > , the result for Qt depends on how the charge 
velocity varies along the wire, and is thus non-universal. On the other hand, the 
charge conductance is universal. Therefore, their ratio is non-universal and the 
Wiedemann-Franz law is violated. 

In a step-like model of Sec. 7.3.1, the transmission coefficient of plasmons is 
equal to 

= — r^FFTi^ — r~- 

J- -t- t'ln 

Obviously, |i (0) |^ = 1 regardless of K, an agreement with what was said above. 
For T <C Wi, the Lorentz number is close to the universal value of 7r^/3e^. For 
T ^ ul, the oscillations of \t {uj)f become very fast, so that |f (w)]^ can be 



replaced by its averaged value 



(I* Ml ) = / —I* Ml 

Jo 



The thermal conductance increases linearly with T, so that Lq approaches a con- 
stant but a non-universal value 

2K 

L\t»^l = -j^2^^o < ^0- (7.27) 

As the Lorentz number varies with temperature in between two limits specified 
by Eqs.(7.26) and (7.27), the Wiedemann-Franz law is violated. 
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Appendix A. Polarization bubble for small q in arbitrary dimensionality 

The polarization bubble in Matsubara frequencies and at T = is given by 

n(«w,g) = ^.^D+i J J d°pdeG{ie + iuj,p + q)Giie,p) 

^ ^ d^pde- — [G {is + iLJ,p + q) — G {ie,p)] 



Ns f ^Dpfi\P + Q])-fip) 



n{itj,q) = -NsVd j ■ 



where / is the Fermi function. Expanding in q, and switching from the integra- 
tion over d^p to d^, we obtain 

(Kl iui 
Q.D ioj — vpqcosO ^ 

where ilo = 47r (3D), = 27r {2D) , = 2 (ID) and ud is the DoS in D dimensions 
per one of the Ng isospin components. For D= 1 , the integral over ft is understood 
as a sum of terms with cos 6 = ±1. It is obvious already this form that the small 
g— form of the bubble depends on the combination ui/vpq for any D. The final 
result depends on the dimensionality. Performing analytic continuation to real 
frequencies iuj ^ oj + iO"*", we obtain 



\ J Ud CO — vrqcosO + tO J 

Taking the imaginary part 

——S{uj — VFqcos6). (A. 1) 

From here 

cos^ = u/vpq, 

which means that 9 w 7r/2 for oj <C vpq. Thus, the fermionic momentum^ is 
almost perpendicular to the bosonic one, q, in this limit. 
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Appendix B. Polarization bubble in ID 



Appendix B.l. Small q 

Free time-ordered (causal) Green's function in ID is equal to 

1 



Gl {s,k) = 



where 

and ± signs correspond to right/left moving fermions. We will be measuring the 
momenta from the corresponding Fermi points. For +branch -.k — kp ^ k and 
for -branch: k + kp ^ k. Consequently, 

Glis,k) 



e^VFk + iO+sgnk 



I assume the Ng— fold degeneracy (Ng = 2 for electrons with spin, A^^ = 1 for 
spinless electrons), so that 

n± (w, q) = '^—Ns [ de I dkGl {e + Lo,k + q) G^ (e, k) . 

(27r) J J 

Calculate, e.g., n_|_ : 

i f , f „ 1 1 



U+{Lj,q) = --l^Nsjdejdk 



(27r)^ J J £ -\- uj — VF{k + q) + iQ'^s,gn{k + q) £ — vpk + iQ'^sg'ok 

= I d£ I dk \- ^ (B. 1) 

(27r) J J to - vpq + iO+sgn(fc + qj - iO+sgnk 

X [Gl (e, k)-Gl{e + uj,k + q)] . (B. 2) 

The integral of the Green's over the frequency gives a Fermi distribution function 
[23] 



n+{k) = -i J ^Gl {£,k). 



For free fermions, 

n+ {k)=e{-k) 

Now, 



no (w, g) = ^ / dk — — — ^ — — — - \e i-k) -0{-k- q)] . 

+ ^ ' 2n J uj-VFq + iO+sgn(fc + q) - iO+sgnfc ^ ' ^ '^^J 
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The integral is not equal to zero only if the arguments of the 6 -functions are of 
the opposite signs. Consider different situations. 

1) k > Q;k + q <Q ^ Q< k < -q ^ q <Q.ln this case, 

no (a;,9) = ^ 5<0; 

2) k<Q,k + q>Q^-q<k<Q^q>Q 

Combining the results for g > and q <Q together, 

n^(^,9) = ^ h?^ — • 

^ ZTT w — Uirg + zO+sgng 



Similarly, 



The total bubble 



n° {u:,q) = ^— . (B. 4) 

ZTT oj — Vpq + ^0+sgna; 



2 



no (a;, g) = no (a;, 9) + n^ (a;, 9) = , ;\ , ■n+ - 5) 



In what follows, we wiU also need the retarded and advanced form of the 
bubble. These forms can easily be obtained by repeating the calculation above in 
Matsubara frequencies and analytically continuing ^ uo+iQ ■ Even simpler, 
one can use the general relation between time-ordered and retarded propagators 
[23] (which works equally well for fermionic and bosonic quantities) 

n^(w,q') = n± (w, g) , for w > 
= n^ (w, q) , for w < 0. 

Using Eqs. (B. 3) and (B. 4) we obtain 

n^{co,q) = ±^ ^— - (B.6) 

ZTT w — vpq + «0+ 



and 

n« {w,q) 



Ng vpq 



2 



TT Lo"^ — Vpq'^ + iO+sgnco 



(B.7) 
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Appendix B.2. q near 2kF 

We will also need the 2kF bubble. This time, I choose to do the calculation in 
Matsubara frequencies: 

n2feF q) = j '^^ J '^^'^^ + iLO,k + q) G_ (ie, k) 

N. 



dk / ds 



' s + oj + i {k + q) e — ik 



Poles in ei = ik and £2 = —i{k + q)—u) have to be on different sides of the real 
axis, otherwise the integral is equal to zero. Choose q > 0. Then this condition 
is satisfied in two intervals of A; : A; > and —A/2 < fc < —q, where A is the 
ultraviolet cut-off 



n. 



2kF — 



An 



A/2 



In- 



dk 



iAvi 



dk 



A/2 



1 



uj + 2ivFk + ivpq 



ui + ivpq 



In 



(J — ivpq 



—iAvF 



47r + Vpq 



2 " 



(B.8) 



Because the result depends on q^ there is no need for a separate calculation for 
the case q <0. 



Appendix C. Some details of bosonization procedure 

Appendix C.l. Anomalous commutators 

p 

P± = X] 4-9/2^+9/2- 

p>0(p<0) 



The operators of full density commute. The operators of left-right densities have 
non-trivial commutators. For example, let us calculate [p+ (g) , p+ {q')] 
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C++{q,q') = [p+ (q) , p+ {q')] = ^ [o'l-q/2%+q/2,"'l-g' 



/2'^fc+g'/2 



p>0,fe>0 



/ a: 



E 

p>0,fe>0 



'p-q/2 



O-k+q' /2 



—a 



-'5p+9/2,fe-,'/2-"I_<,'/2%+(,/2 



V 



-*fc + g'/2,p-,/2 "p-5/2"je+<,'/2 



The first 5— function means that k = p + q/2 + g'/2 > and the second one 
that A; = p - q/2 - /2. 

C++ {q, q') = «p-g/2ap+9/2+9'^^ (P + 9/2 + 972) - «i_,/2-,'S+?/2^ {p - q/2 - q' /2) 
-[f{q,q')-f{q',q)], 

where 

f{q,q') = E «I-q/24-q'/2«P+9/2«fc+g'/2 

p,fc>0 

p,fc>0 

It is easy to show that / (g, q') = f {q' , q) . Indeed, 
fiq',q) = E «p4«p-9'«fe+9 

p,fe>0 

= re — labelling k < — > p = ^ a^a],aj._g/ap+g 

p,fc>0 

= anticomniuting = ^ a^a^ap+^afc+g, = / {q, q') . 



p,k>0 



Thus 



C++ {q, q') = aJ,_g/2S+9/2+9'^ (P + 9/2 + 972)-ap_g/2-,'ap+9/26' (P - 9/2 - 972) • 

p>0 



Introduce a new momentum 



Q = 



q + q' 
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In the first sum, shift p + /2 — > p and in the second sum shift p — q' /2 — > p. 
Then 

C++{q,2Q-q) = ^ a+.^a^+Q (p + g/2) - ^ (p - 9/2)] 

p>0 

p>-q/2 P>q/-2. 

If the main contribution to the sum is given by the states which Ue either deep 
below or far above the Fermi levels, then the quantum fluctuations in the occu- 
pancy of these states are small, and the operators a}^_Qap+Q can be replaced by 

their expectation values (a]r-Q^p+Q) ~ ^Q,onp = Sq^O {pp — p)- Doing this, 
we find 



(PF PF \ iff 

E - E = ^Qfi2^ 
p>-q/2 p>g/2/ V"' 



Pf rPF \ 

dp — dp 

q/2 Jq/2 J 



qL 

Therefore, 



[P+ {<l),P+{-ci)] = |^,spinless. (C. 1) 



The same procedure for fermions with spin gives 



Similarly, 



[p+,<T {q),p+,a'i-q)] = with spin. 



[P-,<T {q),P-,a'{-q)] = -5<^<7'^, with spin. 



and 

[P+,a {q),p-,ai-'l)] = 0- 

Combining these results together 

[Pa.cr (q) ,Pa',a' (-q)] = aSa,a' Scr,a' T;- , 

Ztt 

where a = ± is the chirality index. For full charge density and current, it means 
that 

27r 27r 27r 27r 
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Similarly, 

[/(9),i^(-a)] = o, 



whereas 



r c / N c / M qV qV qV qV 2 ^ 
[P^ (a), 2;^+ 2^+ 2^+ 2^^ = -.^^. 



In 4-notations, 



TT 



where e^o = e^^ = 0, = -e^o = 1. 



Appendix C.2. Bosonic operators 

Let's check that the representation of density operators via standard bosonic op- 
erators does reproduce commutation relation for density. Expand the density 
operators over the normal modes 

q<0 

where, without a loss of generality, Ag can be chosen real and even function of 
q. Fourier transforming p+ (a;) 

P+{q) = r dxe-^^-jY,Ag.{hg.^^'- + b\,e-^^'A, (C. 2) 

q'>0 

= Ag(9{q)bg + e{-q)bl^y 
Choose q> and substitute (C. 2) into the conomutation relation 

[p+{q),p+{-q)] = Al[bg,bl]=Al = ^^^ 



qL 
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Appendix C.2.1. Commutation relations for bosonic fields and i? 
Using 

-cx><g<cx) V^m^ 
-<x)<9<cx> V^l«l-^ 

we find 



[^{x),^'ix)] = *E^ 



Appendix C.3. Problem with backscattering 

As it was pointed out in the main text, straightforward bosonization of the Hamil- 
tonian for the spinless case encounters a problem if one tries to account for 
backscattering. As backscattering (gi) is just an exchange process to forward 
scattering of fermions of opposite chiralities (32), the Luttinger liquid parameters 
with 51 7^ should be obtained from those with gi = by a simple replacement: 
92 — > 52 — 51 • However, if we do this, we cannot satisfy the Pauh principle which 
says that for a contact interaction, when = 54 = 51, all the interaction effects 
should disappear. Indeed, Eqs.6.2) and (6.11) for u and K, correspondingly, 
change to 

2 /1 , 34^2 (92- gi^"^ 



U 



fl + ^ 



27r 



94-92+91 
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For contact interaction, when g2 = gi^ = g\, we get 

2 



K = 1. 

The charge velocity is different from 1. In addition, the product uK is renormal- 
ized from unity-this is also a problem, as it means that the current operator is 
renormalized by the interactions. How to fix this problem? Ref. [77] shows how 
to arrive at the expressions for u and K which satisfy all necessary constraints 
just on the basis on Galilean invariance and dimensional analysis. Ref. [78] ar- 
rives at the same result by using a careful point- splitting of the operators. Here, I 
present the method of Ref. [78]. 

Recall that the density operator, represented in terms of bosonic fields, con- 
tains not only the lowest harmonic (q 0), corresponding to long-wavelength 
excitations, but also harmonics oscillating al q = 2kF , Akp , etc. Indeed, taking 
into account only the 2fcir— oscillations, we have 

p (x) = (ijl (x) e-''""' + (x) e*'^^^') (V'+ {x) e''""' + {x) e'^'^^^) 

= V+ (2;) 1/-+ {x) + -0^ [x) [x) + e-2*'=^^7/'+ (2;) ip- [x) + H.c. 

The first term in this equation has to be treated using the point- splitting proce- 
dure, because it involves two fermionic operators at the same point. The result is 
an infinite constant, po , which is just a uniform density, plus the gradient term. 
The 2kF -component can be bosonized without a problem, as it involves products 
of different femnions. The result is 

p ix) - po = -^dx^p + -J— exp \2y/TTip + 2kFx\ -\- H.c. 
Using this expression for the interaction part of H, we have 

-f^int = \ j '^^ j ^^'^ ~ IP (^) ~ Po\ [P i^') - Po] 

= Hp + Hb , 

where the forward and backscattering parts of the Hamiltonian are given by 

Hf = 2~ y '^'^ y '^^'^ ^ •'^') '^x'-pdx'f, 
Hb = l[^)^JdxJdx'Vix-x') 

X {exp [2iV^ip (x)] exp [-2iV^ip {x')] e^^fe^C^^-^') + H.cj . 
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In Hb, we neglected the terms that oscillate with x, a;',and x + x', and kept only 
those terms that oscillate with x — x'. As our potential is sufficiently short-ranged, 
the oscillations of the first group of terms will average out, whereas the second 
group will survive. Introducing new coordinates 

_ x + x\ 

R = 

r = X — x , 

and assuming that |-R| ^ \r\ , the forward-scattering part of the Hamiltonian 
reduces to 

HF = ^JdR (dRipf J drV (r) = ^ j dR {dn^f . 

The product of the two exponentials needs to be evaluated with care. Applying 
the Baker-Hausdorff identity 

we get 

exp [2i^/TT(p (x)] exp [— 2i Vttv' i^')] = sxp [2i^/n {(f (x) — (p (a;'))] : 

X exp[47r(<p {x - x') ip (0) - ip^ (0)], 

Using the expression for the free bosonic propagator 

l^^ix-x')^{Q)-^'m = ^ln^^^, 

and expanding in r = x — x' under the normal-ordering sign, we obtain 

1 

exp \2i^fKip {x)\ exp \—2i\fKip (x')] = — -47r (dmp) r^—2= "^tt (dRip) a^. 

(While expanding, we neglected the first derivative of ip which can be always 
eliminated by choosing appropriate boundary condition.). Hb reduces to 



Hb = "2 (^^j 2na^ j dRidRip) j drV {r)2cos2kFr 
= j dR{dRipf j drV {r) cos2kFr 
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Therefore, the bosonized form of the total Hamiltonian 

V{Qi)-V{2kF) f ,.,2 



Hint — 



1-K 



J dRidRifif 



manifestly obeys the PauU principle. The Luttinger-Uquid parameters are now 
given by 
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